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DETECTION ELECTRONS THE VERY LOW 
ENERGY 


ABSTRACT 
détection des électrons trés basse énergie dans spectrométre 

est réalisée par d’un accélérateur compteur Geiger fenétre 

trés mince. L’ensemble permet détecter avec une transmission 50% des 

électrons d’énergie pratiquement nulle, cela sans les impulsions parasites 

habituelles dans tels 

INTRODUCTION 

The window which prevents the outflow gas from counter 
absorbs low energy electrons. improve the yield one may use very thin 
window give the electrons additional energy before they hit the counter. 
Both techniques are commonly used the nuclear field and have 
been often discussed (Cook 1954; 1953). However, even without 
window, the flowing gas absorbs most electrons energy less than kev. 
the other hand, the rise parasitic background when accelerating 
device used limits measurements energies the order kev. 
Most workers electron spectrometry take advantage the properties 
lens spectrometer study the low energy field. should noted that the 
detector then point counter; the acceleration the exit slit takes place 
between two circular apertures, and large solid angle offered for detecting 
the secondaries that may produced the exit ring. Moreover high 
acceleration (around kev.) preferred with relatively thick window 
(300 avoid the necessity controlling large gas leakage. Such 
large voltage may favor field emission (Hayward 1953) production soft 
X-rays (Cook 1954). 
detector using acceleration for 180° nuclear spectrometer has been built. 
Its characteristics and behavior will now described. 


APPARATUS 


The nuclear spectrometer which the 180° homogeneous field type has 
rectangular exit slit cm.); the detector counter 


received October 1955. 
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Establishment, Valcartier, Quebec. 


153 


¢ 
| 

a 

ay 

al 


154 CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


with rectangular window cm.) few millimeters behind the slit. 
this geometry, the counter window subtends small angle for the drifting 
secondaries. The acceleration voltage low (less than kv.) with thin 
window thick. 

The counter (Fig. has the following dimensions: diameter, 
0.6 cm.; length, 6.0 cm. The extremities are closed hypodermic needle 


INSULATING 


BOX(LUCITE) GLASS 


EXIT SLIT COUNTER SLEEVE 
BAKELITE 

COUNTER 
EXTREMITIES 
O-RING 
ACCELERATING SCALE 
VOLTAGE COUNTER 
DOOR 


Fic. counter with the Lucite insulating box, counter extremities, and 
counter door with the half O-ring. 


glass sleeve. This sleeve made stretched capillary sealed the needle 
with microflame. The tungsten wire has diameter 0.003 in. The total gas 
pressure cm. (argon, 4.8 cm.; alcohol, 1.2 cm.). The counter has 
plateau with slope 0.1% per volt. The background counts per 
minute. The entrance slit along the counter covered stack three 
layers Zapon Aquanite films (Backus 1945) picked from cold water surface 
with flat rectangular rack. half sectioned O-ring protects the window and 
assures good vacuum seal. 

Lucite box (minimum thickness 1/16 in.) insulates the counter from the 
spectrometer and holds the exit slit the spectrometer. The alignment thus 
easily done and the detector itself the accelerator. 

The accelerating potential obtained from high voltage supply variable 
from kv. The voltage across the counter electrodes must constant 
over long period because the slope the plateau. This was made possible 
using voltage stabilizer (Victoreen VXR-1500) across voltage divider. 
Any voltage between 300 and 1000 may applied and kept stable within 
for any length time any change the accelerating potential. Com- 
parisons between readings with and without accelerations are thus facilitated. 
The signal from the counter brought toa preamplifier through well-insulated 
capacitor. 

drift the plateau induced variations the gas pressure. Pressure 
stability was obtained with Cartesian manostat (Gilmont 1951) whose 
float mechanically closes opens electrical circuit. leak then opened 
between two large reservoirs, the first connected the counter and the second 
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containing the gas mixture higher pressure. The pressure kept stable 
within 0.1 mm. checked against differential manometer. 


CALIBRATION 

From the measured intensities without acceleration and with 
acceleration possible determine the transmission the counter 
window function the energy the electrons. The absorption the 
unaccelerated electrons begins given emission energy But since the 
emission energy interval the accelerated electrons are not absorbed, 
the ratio the intensities, with and without acceleration, gives directly the 
transmission the window. the next interval emission energy, from 
even the accelerated electrons are partially absorbed; but 
from transmission measurements the preceding interval, possible 
correct the accelerated electron intensity its real value and then determine 
the transmission the window this second interval from the ratio the 
unaccelerated electron intensity the corrected accelerated electron intensity. 
For other energy intervals, the procedure the same. our case, the source 
was the active deposit thoron and the accelerating voltage was 4.5 kv. The 
correction made for the partial absorption the accelerated electrons only 
for 200 gauss-cm. Fig. shows the intensity curves for Zapon 
window thick. 


200 


100 


Fic. Transmission the counter window. Curve the counting counts per min.) 
for the accelerated electrons against their values. Curve the counting for the unacceler- 
ated beam. The transmission the ratio 
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The same transmission shown Fig. with the abscissa energy. 
shift the left the amount acceleration gives the actual yield the 


TRANSMISSION 
100% 


ENERGY 


Fic. Yield with acceleration. Curve the transmission curve from Fig. with the ab- 
cissa energy units (kev.). Curve gives the transmission for kv. 


system. For acceleration 5.0 kv., 50% for electrons energy close 
and 90% for those 2.5 kev. energy. 
BACKGROUND WITH ACCELERATION 
The cutoff absorption the window around 2.5 kev. Fig. shows the vari- 


ation the background with acceleration. The point rise the curve does 
not seem depend the rate the gas leakage does the intensity the 
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Fic. Background count (per min.) with the acceleration (A). 


count. grid was interposed between the exit slit and the counter. Fig. 
graph the variation the background against the grid voltage for 
acceleration kv. 

possible explanation these curves follows. Secondaries produced 


ann 
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4000 


Fic. Variation the background against the grid voltage (v.). 


the grounded surroundings are not detected since possible choose 
accelerating potential well beyond the cutoff the window, yet not region 
where the parasitic background significant. These electrons however can 
ionize the gas issuing from the counter. The new electrons will penetrate into 
the counter they receive enough energy overcome the total absorption 
limit. 

The density the gas issuing from the counter decreases very rapidly the 
distance from the counter increases. The probability ionization also falls 
very rapidly. The background can thus start given potential rise very 
rapidly with increasing acceleration. 

voltage +1000 the grid leaves the voltage distribution undisturbed 
(Fig. 5). The background low. Any decrease this voltage equivalent 
moving towards the counter the critical distance for detecting electrons from 
ionized particles. The background thus increases toa point where secondaries 
from the grounded slit are all repelled the grid v.). There there 
sharp break the curve. For lower grid voltage, the grid itself emits second- 
aries. 

The actual apparatus did not permit complete verification these 
hypotheses. 

Different metals exit slit did not seem make much difference the 
background count. interesting note that sudden application kv. 
accelerating voltage followed definite increase the background with 
time. 

special apparatus now being built study the performance the 
accelerator and its optimum operating conditions. These results are not yet 
available. 


1000 
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CONCLUSION 


final evidence the ability this equipment detect very low energy 
clectrons the radioactive source was replaced emitting filament biased 
negatively with respect slitted plate. The distribution these electrons 
was then analyzed with the spectrometer. Curves this type were obtained 
with the bias small This has permitted direct verification the 
calibration the spectrometer the lower energy range. 

The detector has also been used for measuring the low energy distribution 
the active deposit thorium (Fig. 6). The detector-accelerator can used 
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Fic. Low energy electron distribution the active deposit thorium. 


for the complete spectrum the Mev. region with minimum corrections 
from kev. This correction easily obtained from the transmission curve. 

The detector thus valuable tool for studying the complete nuclear 
spectrum down the lower energies. Its properties suggest review the 
backscattering and self-absorption characteristics the sources this region. 
may also used for measuring directly low energy spectrum varying 
the absorption edge with different accelerations. 
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ATTENUATION MC. SOUND ALUMINUM 
LOW 


ABSTRACT 


The attenuation sound frequency megacycles has been found 
reach maximum aluminum 155°K. Combination this result with the work 
Bordoni, who obtained maximum kilocycles and 100°K., gives 
activation energy agreement with the value calculated Mason for dis- 
location relaxation process. The maximum greatly increased plastic defor- 
mation. The velocities sound for annealed and cold-worked aluminum show 
maximum difference about 155°K., again indicating dislocation mechanism. 
The maxima are superimposed background attenuation mostly caused 
scatter grain boundaries and diffraction. 


INTRODUCTION 


Much interest has been shown recently the results measurements 
Bordoni (1954) the attenuation and velocity sound metals over the 
temperature range 4°K. 300°K., the frequencies used being 
The velocities lead, copper, aluminum, and silver were found increase 
the temperature decreased, the velocity changes vanishing high power 
temperature near the absolute zero. The attenuation falls exponentially with 
temperature from 300°K. 200°K. and then rises maximum tempera- 
ture characteristic the metal. Below the peak temperature the attenuation 
again decreases with temperature approaching the value zero high power 
temperature near the absolute zero. 

Mason (1955) has demonstrated the relaxation nature the peaks 
attenuation, and combining results for lead with the high fre- 
quency single crystal measurements Bommel (1954) has deduced acti- 
vation energy for the process. Using model considered first Koehler 
(1953), Mason showed that the peaks may well due the displacement 
dislocation from one atomic line adjacent one against the limiting shear 
stress the crystal. From the model determined the limiting shear stress 
the four metals used Bordoni. 

The present work describes attenuation and velocity measurements made 
polycrystalline aluminum over the temperature range 4°K. 
attenuation maximum found 155°K. The magnitude the peak 
dependent the previous thermal and mechanical treatment the specimen. 
Both the velocity and attenuation measurements tend confirm the dislo- 
cation relaxation nature the maximum. 


EXPERIMENTAL METHOD 


pulse technique was used the present experiments. The pulse generator 
produced Mc. pulse several microseconds duration, which was impressed 
X-cut quartz crystal for longitudinal waves, contact having been made 
between the crystal and the metal specimen means Dow Corning Fluid 


received September 26, 1955. 
Contribution from Royal Military College, Kingston, Ontario, Canada. 


159 


160 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 


D.C. 200 viscosity 12,500 centistokes 25°C. The same crystal acted 
receiver the pulses reflected from the opposite face the specimen. After 
amplification the reflected pulses were displayed oscilloscope whose 
sweep was synchronized with the pulsing oscillator. Fig. shows the block 
diagram the circuit. 


PULSE 
CIRCUIT 


AMPLIFIER 


SYNCHRONIZER 


SPECIMEN 
SWEEP 
CIRCUIT 


Fic. Circuit for the measurement velocity and attenuation. 


The coupling fluid, when cooled slowly through its solidifying temperature 
(about showed tendency lose cohesion with the polished face the 
aluminum sample. This did not happen rapid cooling. Polishing the metal 
finer than 3/0 emery gave good contact all temperatures. 

The quartz was unsilvered and shaped electrode, 
recommended Mason and McSkimmin (1951), employed excite the funda- 
mental mode. 

The amplitude successive reflections from the under side the metal 
specimen was found fall exponentially. and are the amplitudes 
the and reflections, then: 


composed two terms; due the attenuation per centimeter path 
through the solid and due seal losses and reflection. Thus 2C/+D 
where the total length travelled pulse through the specimen. 

the usually defined decrement (energy loss per cycle) /2(stored energy), 
then where the frequency and the speed the sound wave. 
For comparison with internal friction measurements made other ways the 
correspondence between the present measurements and values is: 


(at room temperature where 6.3X105 cm./sec.). The low temperature 
measurements were made with the crystal and specimen supported Collins 
helium cryostat. Temperatures were determined using 
thermocouple soldered the aluminum specimen. The temperature readings 
are therefore somewhat inaccurate extremely low temperatures. 
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Relative measurements only velocity were made. This was done 
photographing the oscilloscope trace and comparing the distances the film 
between five pulses. 

The amplitude oscillation could varied over about 50% changing 
the voltage the condenser which discharges through tuned circuit con- 
taining the quartz crystal. This does not change the shape the pulse the 
crystal, which would course result from simple resistance attenuation the 
Output stage. 

The measurements reported were all made the same sample aluminum 
purity 99.99,%. The grain size the metal was about mm. 

EXPERIMENTAL RESULTS 
Figure shows the variation attenuation decrement with temperature 


over the range 250°K. 4°K. The three graphs were obtained the order 
The aluminum was annealed for hr. 520°C. and then run 


DECREMENT 


0.36 


TEMPERATURE 


Fic. Attenuation sound aluminum. 
Cold worked 1.4%. 
obtained. The specimen was then compressed beyond its elastic limit, the 
permanent strain being 1.4%. Run was then made. recovery anneal 
180°C. for hr. was then given and graph measured. change grain 
size was observed. 
Three features the graphs are apparent: 
(a) the background attenuation very high and increases about 0.04 
the temperature decreased, the attenuation approaching constant 


value very low temperatures; 


0.48 
0.36 
@ 
0.36 
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attenuation reached temperature 155°K. for all three 
conditions; 

(c) the magnitude the increase 155°K. markedly dependent the 
previous history the specimen. Cold working (curve has increased 
the attenuation peak over the annealed (A) condition, while the low 
temperature anneal (C), which followed, has reduced the magnitude 
the peak value intermediate between and 


Figure shows the variation velocity with temperature for metal 
specimen conditions and The measurements are relative only. Since 
the work Bordoni has indicated that the relaxation times for dislocation 


>2 
all 

5 

fe) 50 100 150 200. 250°K. 
TEMPERATURE 


Fic. Temperature dependence the velocity sound aluminum. 

Annealed hr. 520°C. 

Cold worked 1.4%. 
movement are much too long very low temperatures for the dislocations 
moved megacycle frequencies, thus the elastic moduli for the cold- 
worked and annealed specimens should identical those frequencies and 
temperatures. have therefore brought the velocity values for the annealed 
and cold-worked specimens together 0°K. 

The velocity sound the cold-worked specimen noticeably below that 
the annealed specimen all temperatures about 250°K. may also 
noted that the two curves are farthest apart temperature approxi- 
mately equal that which the attenuation maximum takes place, i.e. 
155°K. 

The velocity approaches constant value 0°K. This point will not 
pursued but consequence the Third Law Thermodynamics. 

The amplitude variation attenuation was determined certain tem- 
peratures over the complete range, and was found zero within the 
experimental limits. The strains applied were less than 


DISCUSSION 


Background Attenuation 
(1) Attenuation Due Thermoelastic Damping (Zener 1948) 


This arises because the fact that unidirectional stress applied 
polycrystalline solid will encounter crystallites different orientation and 
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thus strain these crystals different amounts. Thermal currents then flow 
across the grain boundaries and the period the applied stress falls between 
the thermal and adiabatic limits defined the grain size and thermal con- 
ductivity then irreversible conversion mechanical thermal energy will 
take place. This will appear increased attenuation the sound wave. 
The thermoelastic damping decreases with fall temperature and vanishes 
near the absolute zero. Bordoni has found evidence this effect 
between 250°K. and 4°K. However the magnitude the effect 250°K. was 
only about decrement i.e. about attenuation Thus 
thermoelastic damping, while certainly present, contributes very 
fraction the attenuation experimentally found. 


(2) Damping Due Viscous Slip Grain Boundaries 


(1947) has found evidence attenuation low frequencies due slip 
aluminum grain boundaries. The effect should become importance above 
250°K. and probably accounts for the suggestion increase attenuation 
the highest temperatures. The smallness the effect lower temperatures 
would preclude its recognition this present work. 


(3) Damping Due Scatter Grain Boundaries 


This has been studied aluminum Mason and McSkimmin (1951) and 
Roth (1948). arises because the anisotropy mechanical properties 
this metal, that sound pulse encounters change elasticity every 
grain boundary. This gives rise scattering, the frequency dependence 
attenuation being apparently governed the relative sizes crystallite and 
wavelength. For grain sizes about one-third the wavelength Mason and 
McSkimmin find that the attenuation varies while Roth finds little 
attenuation dependence for grain size above the present case 
the grain size about The scattering, believe, should then show little 
dependence wavelength and from Roth’s work about db./ft., 
0.2 


(4) Attenuation Due Diffraction and Other Losses 

Since our main interest this work lies the measurements the damping 
maxima, attempt was made this stage evaluate the losses due dif- 
fraction, alignment, and reflection the seal. The diffraction loss probably 
quite high (~0.35 B), the ratio crystal diameter wavelength being 
only 10. Since the wavelength increases the temperature lowered and then 
flattens off toa constant value, the loss diffraction should behave similar 
manner. should noted that the measured increases constant value 
lowered. Reflection losses the seal are small (~0.01 com- 
parison with the background attenuation while alignment losses our 
machined specimen should negligible. Reflection losses the under side 
the specimen which contact with helium gas are also negligible. 


thus suggested that the major causes the background attenuation 
are scatter grain boundaries and diffraction, the variation with temperature 
being chiefly due diffraction loss changes arising because the change 
sound velocity the metal. 
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Attenuation Maxima 155°K. 

Such maxima, reported Bordoni (1954), have been shown Mason 
(1955) due dislocation relaxation process. Dislocations are supposed 
pinned down irregular intervals impurity atoms giving loops average 
length Thermal agitation causes the loops displaced adjacent 
minimum energy positions against the limiting shear stress and the force 
stretch dislocations. Mason calculates that the energy changes associated with 
these two restrictions are respectively: 


where the atomic spacing along the glide plane, the average loop length, 
the shear elastic constant, and the limiting shear stress. The decrement 
resulting from this relaxation then given reaction rate theory 


where the total number loops per cc., and the 
fraction the total loop covered one kink bent portion the dislocation 
connecting the two minimum energy positions. single sized loop would then 
result attenuation vs. temperature curve showing maximum 
specific temperature. From and Bommel’s results Mason computed 
the activation energy for lead. Using this value and the expressions 
for and deduced the values for the four metals lead, aluminum, 
silver, copper. 

higher temperatures the dislocations become freed from the impurity 
atoms and the attenuation rises exponentially. The results shown Figs. 
and tend confirm the suggestion that around 155°K. dislocation movement 
plays important part the strains caused the high frequency stresses. 
Fig. shows that the attenuation depends the number dislocations 
present. Fig. shows that the temperature where the dislocation movement 
greatest the effective elastic moduli have been reduced the cold-worked 
specimen relative the annealed specimen which contains fewer dislocations. 
This shown the fact that the velocities the two specimens differ most 
around 155°K. 

assume our work that the peak damping above the background 
arises from single simple relaxation process and that and are the two 
temperatures which the damping half the peak value, then the activa- 
tion energy given (Zener 1948) by: 


From Fig. curve and are respectively 182°K. and 128°K. whence 
2300 cal./mole which good agreement with Mason’s calculated value 
2260 cal./mole. 

Bordoni observed peak 100°K. for aluminum using frequency 
can also deduce the activation energy from his result and from 
our own with peak 155°K. and frequency 
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where and are the peak temperatures frequencies and This gives 
2700 cal./mole which reasonably close the above results taking into 
account the difference the sample impurity and condition and the difficulty 
spotting exactly the peak temperature. 

Thus the results give strong confirmation the correctness the steps 
theory leading from the Bordoni observed lead results the deduced 
aluminum values. 


CONCLUSION 


The attenuation Mc. sound waves polycrystalline aluminum has been 
found reach maximum 155°K. The activation energy for the process 
agrees with that deduced Mason from theory dislocation relaxation. 


This work was supported Defence Research Board Canada, grant 
number 7510-14, Project D44-75-10-14, whose assistance gratefully ack- 
nowledged. 
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ETUDE EXPERIMENTALE DIFFRACTION DES 
MICRO-ONDES PAR DES OUVERTURES SYMETRIE 


RESUME 

Des mesures diffraction dans région Fresnel par des ouvertures 
révolution ont été réalisées une longueur d’onde 1.25 cm. Les 
ouvertures sont trou circulaire, couronne réseau couronnes 
Soret. Les résultats obtenus montrent que calcul diffraction 
sur Kirchhoff donne des résultats précis mieux 
pour des plans mesure aussi prés que deux fois diamétre 
Cette précision pour des plans plus éloignés. 


INTRODUCTION 


nombreuses études sur diffraction des micro-ondes parues ces der- 
années ont porté plus souvent sur détermination, théorique 
des conditions électromagnétiques voisinage 
des écrans leurs discontinuités. Dans domaine 
particulier diffraction, rapport étendu publié récemment par Bouw- 
kamp (1954) fait une mise point détaillée des progrés réalisés des résultats 

Une approche différente est celle plusieurs chercheurs qui 
ont dans les micro-ondes outil travail remarquablement commode, 
rendant possible tout moins plus aisée d’une multitude 
problémes appartenant propre point vue est notamment 
celui Toraldo Francia (1952), Carrara, Checcucci Schaffner (1948) 
Italie, comme aussi celui, ailleurs, Hull (1952), Kastler (1952) Zucker 
(1952). diffraction par exemple, les minuscles systémes franges observés 
dans région optique deviennent avec les micro-ondes trés gros phéno- 
méme des dimensions usuelles. L’existence des moyens 
détection appropriés fait immédiatement penser transposer profitablement 
hyperfréquences photométrie, pénible, des taches diffraction op- 
tique. cette peut espérer obtenir avec les micro-ondes une infor- 
mation détaillée sur des questions grande importance théorique, comme 
détermination précise domaine des approximations admissibles, tant 
celles qui relévent simplement géométrie que celles, plus fondamentales, 
sur rapport (dimension des ouvertures/longueur dont dépend 

présent travail constitue effort dans cette voie. cherché d’abord 
mettre point dispositif expérimental diffraction des micro-ondes 
permettant des mesures fines reproductibles des diagrammes diffraction. 
recu pour premiére fois février 1955, sous forme revisée novembre 
Contribution Département Physique, Faculté des Sciences, Université Laval, Québec, 


Québec. 
Canadian Armament Research and Development Establishment, 


Valcartier, Québec. 
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L’installation une fois point obtenu une série courbes expérimentales 
d’un bon nombre diffraction révolution pour plupart 
desquels des calculs optiques étaient disponibles. Spécifiquement, étudié 
diffraction Fresnel par des ouvertures circulaires, des ouvertures annu- 
laires des réseaux annulaires. 

Dans une premiére partie présent mémoire nous passerons rapidement 
revue les travaux théoriques cause. Une description 
diffraction aux micro-ondes réalisée suivra. Nous donnerons ensuite certain 
nombre des courbes diffraction obtenues nous terminerons par quelques 
conclusions générales qui ressortent notre étude. 

REVUE SOMMAIRE DES THEORIES SOUMISES VERIFICATION 

2.1. Parmi diversité des phénoménes diffraction, ceux qui pré- 
sentent symétrie révolution forment une classe d’un intérét spécial 
raison des applications optiques. Considérons écran révolution éclairé 
par une source ponctuelle radiation électromagnétique placée sur 
observe dans des plans normaux systéme complexe franges cir- 
culaires centrées sur méme axe. L’ensemble réalise ainsi maniére 
‘‘a Si, dans arrangement équivalent, systéme 
comporte une lentille supposée parfaite (ou plusieurs), les franges Fresnel 
peuvent encore observées dehors plan focal, sont alors connues 

L’écran révolution plus simple est une ouverture circulaire pratiquée 
dans écran plan opaque son complémentaire, disque opaque. L’un 
ces obstacles représentent des cas limites d’une ouverture annu- 
laire plan révolution plus général consiste systéme 
d’ouvertures annulaires concentriques dont les rayons bord sont assignés 
d’une fagon quelconque. Dans cas général, diffractée par 
systéme annulaire point situé hors est obtenue 
superposant les amplitudes contribuées point par les couronnes libres 
individuelles. Or, pour quelconque ces couronnes libres, est immé- 
diat que diffractée s’obtient soustrayant con- 
tribuée par cercle bord extérieur celle qui provient cercle bord 
intérieur. Les problémes diffraction proche révolution raménent donc 
tous définitive celui circulaire, dont solution apparait 
ainsi fondamentale. 


2.2. Diffraction par une ouverture circulaire 

Les franges diffraction proche circulaire ont été beaucoup 
étudiées depuis Fresnel Poisson. Dans mémoire classique 
fondamental Lommel (1885) parvint premier calculer répartition 
dans systéme complexe franges. Peu aprés cependant 
Struve (1886) obtenait des résultats semblables quant principe des calculs 
mais moins détaillés. 

figure illustre géométrie probléme. Soit point quelconque 
coordonnées polaires dans plan diffringente circulaire 
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distance correspondante point point d’observation dernier est 
situé sur cercle rayon centré sur D’autre part soit rayon 
l’onde sphérique incidente, distance jusqu’au plan d’observation 
monochromatique utilisée. Par suite symétrie révolution suffit 
considérer les phénoménes dans plan horizontal XOZ. 

réalisée, diffractée est donnée par étendue 
l’ouverture 


a 


Dans cette expression 

étant supposée constante sur front d’onde sphérique incident 
représente les premiers termes développement série chemin optique 
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convient position point d’observation moyen des 
paramétres 


Introduisant ces derniers intégrant par rapport obtient 
a 
0 


intégrale transcendante d’évaluer que nous désignerons par 

2.21. d’un processus d’intégration par parties, Lommel parvint 
représenter fonction moyen deux développements série, fonc- 
tion des paramétres ka? Ja. Ces quantités ont sens physique 
important, premiére représentant double différence phase entre 
chemin optique axial chemin optique marginal. Quant vaut 
moitié différence phase pour deux rayons arrivant point 
provenance des deux extrémités d’un diamétre 

Lommel calculé pour circulaire plusieurs diagrammes 
diffraction, soit pour les définis par 27, 
dans méme mémoire probléme disque opaque. plus, fait des 
observations étendues des franges circulaire vue comparer 
les positions des minima des maxima avec les prédictions théorie, 
obtenant dans chaque accord excellent. Par suite, seule étude photo- 
métrique ces franges qui ait été tentée est due Maréchal (1948) con- 


2.3. Diffraction par une ouverture annulaire 

Parmi les résultats concernant les franges défaut mise point pro- 
duites par une ouverture annulaire largeur finie mentionnons les diagrammes 
diffraction pour annulaire rayon bord intérieur égal 
0.707 rayon extérieur calculés par Boivin. Les isophotes voisinage d’un 
foyer parfait pour une ouverture annulaire dont rayon bord intérieur est 
égal dans cas 0.25 rayon bord extérieur dans autre cas 0.707 
ont été calculées par Linfoot Wolf (1953) par une application directe des 
formules Lommel. Les diagrammes diffraction que peut tirer ces 
réseaux d’isophotes manquent cependant précision. 


2.4. Diffraction par 

Réseaux 

principe, calcul direct des expressions existantes pour cas 
circulaire appliquée successivement aux cercles 
annulaire est possible pour petit nombre couronnes concentriques. Mais 
cette tache devient prohibitive nombre des cercles cause 
est peu élevé des difficultés convergence présentent dés que 
paire variables associées chaque cercle prend des valeurs considérables. 
Afin d’améliorer cette situation, Boivin (1952) publié une méthode calcul 


2 


170 CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


des franges diffraction proche, toujours dans région 
Kirchhoff est valable, dues systéme quelconque d’ouvertures annu- 
laires. 

Parmi les types réseaux annulaires étudiées, certains présentent des 
bandes annulaires trop étroites pour que les franges calculées puissent étre 
vérifiées utilement avec des micro-ondes. Dans présent mémoire seuls 
seront présentés discutés des diagrammes diffraction théoriques ex- 
périmentaux obtenus pour des réseaux Soret petit nombre cercles. 
cercle est rayon premier cercle. 


MONTAGE EXPERIMENTAL 


schéma montage est présenté figure tandis que figure 
est une photographie. klystron 2K33 opérant dans bonnes conditions 
stabilité fréquence puissance est comme générateur d’hyper- 
fréquence. Une onde carrée bonne qualité, variant taux 1000 cycles 


Source 
de 
Tension 


diffringent 


coteur 


Amplifi- 
cateur 


par seconde, module fréquence puissance émise sont 
respectivement moyen d’un d’un bolométre 
disposé dans une ligne couplée ligne principale. adapteur d’impédance 
cornet électromagnétique. Les dimensions dernier ont été choisies pour 
obtenir diagramme rayonnement d’intensité constante sur 
étudiée. Les écrans diffringents sont montés centre d’un écran opaque 
dimensions grandes pour prévenir toute interférence nuisible 
entre faisseau étudié radiation passant autour 
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Fic. Photographie montage. 


diffractée est mesurée moyen sonde auto- 
matiquement; captée est indiquée sur enregistreur graphique. 
Cette sonde est demi-onde; celui-ci permet mesure des champs 
diffractés jusque dans les ouvertures sans introduire d’erreurs appréciables. 
est alimenté par conducteur coaxial rigide s’amenant dans 
direction propagation des ondes afin toute perturbation des mesures. 
une distance quelques dizaines longueurs d’onde sonde, coaxe 
est courbé hors champ. L’énergie captée est transmise par conducteur 
coaxial une transition coaxe guide d’onde; dernier est terminé 
par cristal silicium. L’onde redressée, variant taux 1000 cycles par 
seconde, est amplifiée d’un amplificateur sélectif haut gain accordé 
fréquence modulation source. signal audio obtenu, proportionnel 
puissance absorbée par sonde, est finalement redressé appliqué 
cision nous permis vérifier que l’ensemble montage possédait une 
linéarité excellente. 

Pour permettre des mesures automatiques, sonde est montée sur 
chariot déplagant sur banc coaxial rigide, transi- 
tion coaxe guide d’onde ainsi que support cristal sont aussi montés sur 
chariot. Une vis sans fin actionnée par moteur synchrone entraine 
chariot par d’un systéme d’engrenages vis tangentes. banc 
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d’optique est monté sur deux rails paralléles optique 
syst¢me permettant mesurer les figures diffraction dans des plans 
paralléles sur une distance trois métres partir dernier. 
moyen montage, les champs diffractés peuvent étre étudiés dans région 
proche dans région Fresnel, fagon rapide, d’une 
courbe requérant moyenne une minute. 

Afin d’obtenir une distribution d’intensité constante incidente 
sur les écrans diffringents étudiés, est nécessaire d’utiliser des cornets pos- 
des diagrammes rayonnement trés divergents. Ceci diminue 
puissance utilisée augmente radiation frappant les cloisons 
ment revenant, plusieurs réflexions, traverser interférer avec 
faisceau principal diffracté. Pour éliminer les erreurs mesure produites 
par ces réflexions, région voisine contour été couverte 
centimétriques. Pour obtenir une précision mesure 
quelques pour cent, une grande quantité ces absorbants été 
requise. Des mesures réflexion produite par différents matériels nous ont 
conduit réalisation d’un absorbant trés bon marché qui révélé 
idéal pour notre fin. matériel est laine d’acier commerciale; montée sur 
des panneaux carton, elle fournit moyen simple trés convenable 
rendre appartement sans écho. Ses qualités d’absorption réflexion sont 
comparables aux absorbants micro-ondes commerciaux, fréquence 
utilisée (24,000 Mc./s.). 

Les dimensions maxima des ouvertures utilisées sont limitées par longueur 
diamétre ont été étudiées dans région Fresnel dans région 
prés des ouvertures. Des études ont également été faites sur des ouvertures 
annulaires dont rapport entre rayon interne rayon externe est 
0.707, sur des réseaux cercles ayant pour rayon premier cercle les 
valeurs 7.5, 12.5 cm. 

RESULTATS 

Les courbes expérimentales diffraction dans région Fresnel par 
des ouvertures circulaires longueurs diamétre sont pré- 
sentées figure Ces courbes ont été mesurées dans plan dans 
plan sur les graphiques elles sont comparées aux résultats théoriques 
Lommel. pas notre connaissance photométrique d’en- 
semble correspondant aux courbes diffraction Lommel dans domaine 
visible. L’accord entre les courbes théoriques expérimentales est trés bon 
pour les plans mesure éloignés Car alors caractére vectoriel 
p*) font sentir. Mais pour les mesures prés figure 
diffraction mesurée dans plan magnétique est plus étroite que celle mesurée 
dans plan électrique. Ces différences sont attribuables nature sonde 
qui détecte que composante champ électrique paralléle aux bras 
or, cette composante est différente suivant position point d’ob- 
servation égale distance autour différence d’inclinaison des 
rayons ces points. 
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a2 


0.6 


€ 


Fic. Fresnel par des ouvertures circulaires. théorique 
Lommel, expérimentale plan expérimentale plan 1.25 cm. 


Les positions des maxima des minima des courbes diffraction sont bien 
vérifiées plan calcule que les expressions théoriques 
Lommel sont bien vérifiées mieux environ, pour des plans mesure 
situés une distance égale plus grande que deux fois diamétre 
ture. qui apparait remarquable regard des hypothéses fondamentales 

Les figures représentent les graphiques des courbes diffraction par 
des réseaux cercles centre opaque, trois diamétres différents 
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mesurés dans les plans d’observation Sur ces graphiques 
elles sont comparées aux courbes théoriques calculées par Boivin. Onconstate 
encore que les figures mesurées dans plan magnétique sont plus étroites que 
celles mesurées dans plan électrique. Dans cas plus petit des trois 
réseaux cm.) constate, plus qui est mentionné ci-haut, des 
différences dans position des minima des maxima secondaires ainsi que 
dans leur amplitude; différences qui s’accentuent passant plan 
plan 2.47, dire mesure que approche réseau que les rayons 
deviennent plus inclinés sur des formules théoriques 
ces plans est surtout due grande inclinaison des rayons, dont pas 
tenu compte dans Kirchhoff. peut montrer encore dans 
cas des réseaux que les expressions théoriques sont vérifiées mieux 
pour les plans d’observation situés des distances plus grandes égales 
environ deux fois diamétre réseau. 

variation pouvoir résolutif d’un réseau avec nombre couronnes 
été étudiée est montrée dans série courbes présentée figure 
voit que courbe diffraction d’un réseau rapproche celle d’Airy 
lorsque nombre couronnes croit. résolution diminue fonction 
nombre couronnes; pic central plus fin étant obtenu pour réseau 
d’une seule couronne. 

Les figures diffraction d’un réseau centre ouvert cm. ont 
présentée; elle est comparée courbe diffraction d’un réseau semblable 
centre opaque aussi figure d’Airy. constante que résolution d’un 
réseau centre opaque est plus grande que celle lentille méme dia- 
métre, mais que les maxima secondaires sont plus intenses. Par contre, 
réseau centre ouvert posséde une résolution plus faible que cette derniére 
systéme d’anneaux diffraction moins intense. 

CONCLUSION 


Les théories diffraction, basées sur Kirchhoff, 
permettent calcul diffractée avec une bonne précision con- 
dition que plan soit situé des distances égales 
plus grandes que deux fois son diamétre. des points d’observation plus 
prés, les théories récentes Woonton, Bekefi Neugebauer (1952) doivent 
étre utilisées. diffraction dans région Fresnel par des systémes d’an- 
neaux concentriques peut calculée appliquant les expressions Boivin. 
Pour les réseaux Soret, nous avons observé que résolution est plus grande 
que celle circulaire lorsque centre est opaque, alors qu’elle 
est plus faible lorsque centre est ouvert. Cette résolution est maximum pour 
réseau d’une couronne, pour annulaire. 
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SCATTERING ELECTRON ATOM AND 
REARRANGEMENT 


Ta-You 


ABSTRACT 


For the problem the scattering electron hydrogen atom, the 
complete equivalence shown the treatment Trainor and based 
symmetry considerations and that Altshuler giving direct proof the 
method Mott and Massey. generalization the treatment Altshuler 
scattering processes involving non-identical particles other rearrangement 
collisions made. Also, the usual treatment involving matrix elements for 
transitions between states belonging two different Hamiltonians discussed 
and justified. Finally, the recent treatment Dirac rearrangement collisions 
discussed, and shown that correct formulation the method and 
interpretation the result give exactly the same result the usual treatment. 


INTRODUCTION 


well known that the perturbation treatment rearrangement 
collisions, the type C+D, say, where are atoms, 
molecules, electrons, there the difficulty that the initial and the final 
states are not the eigenstates single zeroth order Hamiltonian (Schiff 
1949). completely convincing arguments seem have been given for using 
non-orthogonal states the transition matrix elements (see footnote 3). 
Recently Dirac (1955) has given treatment take into account the effect 
non-orthogonality. This treatment leads different result from that 
the usual treatment even what corresponds the first approxima- 
The purpose the present work (i) start from the rigorous treat- 
hydrogen atom, which shown equivalent much simpler treatment 
Trainor and (1953) based the property the system 
the two electrons, (ii) extend Altshuler’s treatment other rearrangement 
collisions, (iii) justify the usual perturbation treatment involving non- 
orthogonal states, and (iv) show that Dirac’s recent formulation and result, 
while apparently different from the usual treatment, are fact completely 
identical with the usual ones when properly carried out and interpreted. 

The scattering electron hydrogen atom was first treated 
Oppenheimer (1928) the perturbation method. rigorous but incomplete 
treatment the problem was given Mott and Massey (1933), and the 
completion this treatment has been the subject few recent studies 
(Trainor and 1953; Altshuler 1953; Moses 1953). Both because some 
rather important points are not sufficiently recognized and because the 
possibility extending the treatment the problem other rearrangement 
collisions which the scattering electron hydrogen 

received original form August 31, 1955, and, revised, November 16, 1955. 
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atom perhaps the simplest example, shall begin with brief resume 
the treatments Trainor and (1953) and Altshuler (1953) and 
demonstration that these two treatments are exactly equivalent. 

The Hamiltonian the system is, atomic units, 

(1) 

The initial state consists one electron incident with momentum 
hydrogen atom the state The problem find the probability one 
electron being scattered (represented asymptotically spherical wave) 
and the other electron being hydrogenic state 

The simplest treatment this problem make use the symmetry 
property the system the outset (Trainor and 1953). Let 
any function two electrons. always possible expand the 


One can construct the symmetrized (i.e., symmetric antisymmetric) 
function 


The function 

can expanded again the complete set 

that 


(r1) gn (2). 


(6) 


The above holds without any reference any specific problems. 

For the scattering problem, shall seek symmetrized solution the 
Schrédinger equation 


such that represents the initial state (consisting electron 
incident with momentum and the other electron the state plus the 
scattered waves (one electron outgoing and the other state @,), 


tKnTi 


where the energy the hydrogenic state ¢,. Since sym- 
metrized, satisfies the same condition (8) with and interchanged. 
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The Schrédinger equation (7) with the asymptotic condition (8) can 
rewritten the form the following integral equation 
where the function the operator and the 
parameter 


9 9 2 
(11) 
lr.’ A exp (ik, |) * 


and given (9). From (10), one readily obtains the scattered amplitude 
(8) 


9 


which the primes the variables integration (10) have been dropped. 
Equations (10) and (14) are exact the sense that approximations have 
been made. Since the function not known, one makes the usual 
approximation replacing the integrand (14) the initial wave 
function which when symmetrized 


follows that the scattered amplitudes for the symmetric and antisymmetric 
state (singlet and triplet state respectively) are given 

where 


Fis 


These are the results Oppenheimer. That the last two expressions for 
are equivalent follows from general theorem (Schiff 1949) and can readily 
verified noting that the matrix element with respect the 

The difference between the above treatment and the usual one Mott and 
Massey (1933) can appreciated after brief resume the latter. 
describe the scattering, solution the form (2) sought 
that satisfies the asymptotic condition 
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with given (9). Let the function determined the condition 
the form (5) 


(5A) = f G,, bn 
n 
sought which satisfies the asymptotic condition 
iknr2 
(19) G, (fe) £n (Oe, $2). 


able describe the scattering involving two distinguishable electrons, 
necessary show that the two solutions and are 
equal before one can form (16) for the symmetric and the antisymmetric 
states. That 


not itself obvious; for this means that the function determined 
the asymptotic condition (18) for scattering the same time also 
satisfies the asymptotic condition (19) for scattering. the 
preceding treatment Trainor and Wu, the total (direct and exchange) 
scattering obtained together (10) and (14). But proof (20) has been 
given Moses (1953) and Altshuler (1953).-The latter method 
ticularly interesting since can extended the case involving non- 
identical particles example, the collision negative meson and 
hydrogen atom), and since its equivalence the result (10) Trainor and 
can readily shown. shall describe briefly follows. 

The Schrédinger equation for the unsymmetrized wave function 
(2) subject the asymptotic condition (18) for scattering 
replaced the integral equation 


introducing the function 


(23) 


and making use the hermitian properties the operator with respect 
the Green’s function (12) and (22), can shown that (21) 
can transformed the following equivalent forms 
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2 | | 


seen from (24) (25) that satisfies the asymptotic condition (19) 
and fact making the approximation 


the integrand, (24) gives for the scattered amplitude the last expression 
(17). 
The complete equivalence between this treatment and that Trainor and 
seen interchanging the particles and (25), thus 
9 


e lal iT; — Bs | 


and forming the symmetrized function from (21) and (26) 


The result seen exactly the same Eq. (10). 

Equations (21) and (25) state that the function satisfying the asymp- 
totic condition also satisfies the asymptotic condition the form 
(21), the condition for outgoing waves explicitly exhibited while states 
representing outgoing waves are implicitly contained disguised form 
(mostly hidden the terms vice versa for (25). 


REARRANGEMENT COLLISIONS 


The importance the proof that (21) can transformed into 
the form (24) (25) does not lie much enabling one obtain the 
scattered amplitudes (16) and (17), for these could have been derived 
very simply (10) symmetry arguments. The importance lies rather 
the fact that the proof the transformation from the form (21) into the 
form (24) (25) does not depend the non-distinguishability the two 
electrons; depends only the hermitian property the operators (namely 
the Hamiltonians and with respect their Green’s functions. 
This hermitian property turn depends only certain properties the 
functions, namely, that they not diverge infinity and not 
possess singularities order higher than the simple pole. Thus the equivalence 
the two forms (21) and (24) can proved for any rearrangement collisions 
systems for which the functions possess these properties. For the 
problem the collision between negative meson and hydrogen atom, 
the foregoing result denote the coordinate and momentum 
the meson, and the corresponding quantities the electron, 
have only make the following replacements (21) and (24) (25) 


The scattered amplitudes for the direct and the exchange scattering are given 
(17) one makes the replacements 


writer wishes thank Dr. Singh for making explicit calculations generalizing 
the proof Altshuler. 
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where the wave function for mesic atom. 

Similarly, the mathematically equivalent problem the collision between 
positive ion other than proton (say, deuteron) and hydrogen atom 
can treated the same manner. The scattering the ion and the capture 
the electron from the hydrogen atom the ion correspond the direct 
and exchange scatterings the meson hydrogen atom problem. The removal 
the assumption infinitely heavy nucleus the latter problem only 
introduces some further non-essential details the treatment. 

IV. NON-ORTHOGONALITY STATES REARRANGEMENT COLLISIONS 

the scattering problem treated the so-called stationary state 
method. The scattered amplitudes and (17) are obtained from the 
exact theory means the approximation (15). necessary that the 
same approximation, the same result should obtained from the time- 
dependent perturbation theory which the scattering treated transition 
caused perturbation. From the point view, the matrix 
element (17) for the exchange scattering presents some difficulty that 
needs some discussion. g,, the initial and the final states, 
and respectively, are not the eigenstates the same 
zeroth order Hamiltonian. allow for the non-orthogonality these 
states that Dirac (1955) has recently suggested another treatment re- 
arrangement collisions. shall show, however, that this difficulty only 
apparent, and shall discuss Dirac’s recent suggestion the next section. 

Let the Hamiltonian the system grouped two alternative forms 


such that appropriate for the description the initial, and the 
final state, i.e., wish find the transition probability from eigenstate 
(The ambiguity which these two perturbations used only 
apparent and will mentioned again below.) The Schrédinger equation the 
system 


= 
(29A) (i eae = V; y. 
Let expand the the left-hand side the complete set 


making use the relation 
Hy xn Xn 
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one obtains 


Let expand the the integrand above the complete set 


Substitution (32) into (31) gives exact system equations. Let assume 
that initially the system the state 


first approximation, one may neglect the integral 
compared with 


integrate (31) with these values and obtain the usual manner for the 
transition probability per second the system going from the state any 
one the states the continuum the 


where the number states per unit energy range energy equal 
the initial energy the system. can shown general that, subject 
certain hermiticity conditions, 


then seen that (34) and (35) agree exactly with the last two expressions 
(17). 


remains remark that the ambiguity the value the matrix element 
(33) upon changing additive constant account the non- 
orthogonality the functions and can removed choosing the 
Hamiltonians (28) such way that and approach zero 
the separation between the colliding particles (and the resulting particles) 
becomes infinite (Schiff 1949). 

Finally, interest show that theorem essentially equivalent 
that Altshuler (1953) §II can proved this theory transitions. 
The initial condition such (33), together with the time-dependent 
Schrédinger equation (29), completely determines the subsequent times. 
expanding the left (29A) the set (30), one obtains the 


that the use the expansion (32) (31) leads directly (34). Schiff (1949) uses 
the expansion (30) and obtains (34) assuming the possibility the expansion the initial 
state function terms only those the set that have the same energy 
clear that this latter expansion not possible, nor necessary the above treatment. 


186 CANADIAN JOURNAL PHYSICS. VOL. 34, 1956 


coefficients that lead (34) and the (30), denoted appro- 
priate for describing the exchange rearrangement transitions. If, the other 
hand, one writes the Schrédinger equation (29) the equivalent form 


expands the here the set (32), and neglects 
compared with 


then one gets the coefficients that lead the usual result 


| 2 
(34A) = Jot po dr| Pm 


and the solution obtained this way, denoted yp, appropriate for 
describing the scattering. That given initial condition, such 
(33), uniquely determines the solution the Schrédinger equation for sub- 
sequent times requires that the two solutions and equal, i.e., 


VE Vo. 


This true the exact solutions. the first approximation, the result (34), 
(34A) obtained this method exactly that contained (17) given the 
stationary state 


GENERAL REMARKS THE THEORY SCATTERING 


The result (20), namely, that the same wave function describes either the 
stated more general way. The result (20) and its generalization 
scatterings involving two different particles, negative meson incident 
upon hydrogen atom, say, state that the scattered wave may analyzed 
either one two representations: asymptotically outgoing meson 
and hydrogen atom hydrogenic state (discrete continuous), and 
(b) asymptotically outgoing electron and mesic atom 
state (discrete continuous). The totality states the representation 
(a) contains the totality states the representation (b) and vice versa. 
fact, and are not the only representations possible. One might 
wish know the probability that both the meson and the electron come out 
definite directions and energies. For this probability, and are not 
the appropriate representations and one would analyze the scattered wave ina 
representation (c) which asymptotically both particles are outgoing waves. 
Again, one might wish know the probability that the incident particle 
captured inverse Auger transition (radiationless), the final system 
being one the quasi-stationary states (doubly excited, ‘‘discrete’’ states 


‘The writer wishes thank Dr. Allcock for discussion this point. 
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imbedded the continuum). For this information, one would analyze the 
wave function representation (d) which both particles are 
states. All these show that the choice the representation describing 
scattering problem determined the specific information one wishes 
get, and that one uses one representation time. 

from the point view transitions, the theory states that starting 
with given initial state and given perturbation, the state (at time 
the system definite which can analyzed either the set (the 
eigenstates (30), the set (the eigenstates H,°) (32), 
other sets. The total probability for either all the direct scatterings 
all the exchange scatterings (34) is, course, unity according (30) 


(32), 

These equations mean that all the scatterings together contain also 
all the scatterings and vice versa. Using the representation 
(30) does not mean that only direct scatterings are allowed and exchange 
scatterings are excluded. Similarly for the representation (32). One 
cannot obtain without ambiguity the probabilities some direct and some 
exchange scattering the same time using mixture both representa- 
tions, since each state one representation contains many states the other, 
when the two sets and are not orthogonal. 


VI. RECENT TREATMENT REARRANGEMENT COLLISIONS 


Dirac (1955) recently suggested treatment rearrangement collisions 
allow for the non-orthogonality states. Let the Hamiltonian the system 
written the form (28), and let the wave function written the 
form 
where the initial state and the scattered wave that assumed 
small compared with using the relation 

one obtains 


(36) 

Dirac writes known function, the form 

where 


= E oy, 


(38). putting (38) into (36) and neglecting one solves for with 
the condition for outgoing waves. The probability amplitudes for direct and 
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and respectively, and these are obtained from (38), namely, 


dr, 


where 


this formulation, the direct and the rearrangement scattering amplitudes 
and for various states and are coupled together the system 
equations (39) the non-orthogonality integrals obvious that the 
values and obtained solving the system (39), will depend the 
number states and included the sums remove this 
range for the summations (38). For the problem the scattering meson 
hydrogen atom, for example, plausible prescription would seem 
include the sum all the corresponding hydrogenic states that 
are accessible the given initial total energy the system, and include 
all the corresponding energetically accessible states the 
mesonic atom. both sums, the plane waves representing the other particle 
are allowed extend their complete spectra. Momentum conservation 
will taken care the matrix elements themselves and energy conserva- 
tion will come about solving the equation (36). This choice the range 
(38) would give unique value for and but there the 
question already emphasized Dirac (1955) whether the wave functions 
and included the two sums are linearly independent. answer 
this question difficult since involves the question the rank con- 
tinuous infinite matrix. 

Apart from this difficult question linear independence the wave func- 
tions included the two sums (38), there the difficulty this 
method that the transition probabilities |a,|*, given (39) differ from 
those given the treatments the preceding sections, namely, 


Since the result (41) has been shown follow from the exact treatment 
and also justified according the perturbation theory (§IV), 
the result (39) has reconciled with (41). This paradox can resolved 


writer wishes thank Dr. Eden and Dr. Bhatia for discussions this 
point. 
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one extends both summations (38) their complete spectra and expands 
the the second sum (38) the complete set the that (38) 
becomes 


alternatively, one expands the (38) terms the complete set the 
Xu, that (38) becomes 


(42), (43), the summations over and are over the complete spec- 
trum. seen from (42) and (43) that not inter- 
preted the probability coefficient for the direct transition the state 
for the rearrangement transition the state x,. This restores the agreement 
between (39) and the result (41). 

After this reinterpretation (39), one naturally would simply expand 
(38) either the complete set 


(44) Vip = 
the complete set 


putting (44) into the first form (37) one obtains which satisfies the 
asymptotic condition for direct scattering, with the scattered amplitude 
given (41). putting (45) into the second form (37), one obtains 
which satisfies the asymptotic condition for rearrangement scattering with the 
scattered amplitude given (41). The question whether and 
obtained are the same function not answered the affirmative according 
§§II, III, IV. Since contains the outgoing waves the exchange scatter- 
ing implicitly (and contains the outgoing waves the direct scattering 
implicitly), there need for the (38) Dirac. 
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PATTERNS STUB ANTENNAS 
CYLINDRICAL 


ABSTRACT 
Theoretical radiation patterns are presented for radial electric dipoles 
located isolated circular cylinder, anda cylindrically tipped half-plane. These 
cylindrical structures are considered infinite length and perfect 
conductivity. The computed patterns for the isolated cylinder compare well with 
the experimental results reported Bain. 
INTRODUCTION 
The radiation pattern microwave antenna aircraft very much 
dependent the particular shape the surface which mounted 
(Granger and Bolljahn 1955; Wong 1954). Owing great mathematical 
difficulties, exact theoretical analysis not possible for actual shapes. 
For this reason, recourse usually made experiment actual aircraft 
investigations means models (Sinclair, Jordan, and Vaughn 1947). 
The other alternative investigate theoretically the case where the body 
represented some simple geometric shape which amenable rigorous 
mathematical treatment. For this purpose the circular conducting cylinder 
infinite length often chosen with the source antenna being represented 
electric dipole located just outside the 1943; Lucke 1951). 

Unfortunately, the usual mathematical treatment cylinder problems 
electromagnetic theory leads radial mode harmonic series solutions 
which become poorly convergent when the diameter the cylinder greater 
than about one wavelength. possible transform these harmonic series 
solutions, function-theoretic means, so-called Watson residue series 
(Bremmer 1949), which highly convergent deep the region 
for very large cylinders. The series, however, becomes very poorly convergent 
the region near the optical boundary. Furthermore, usual 
geometric-optical approximations also fail give reliable results this 
critical region. therefore seems advisable use the harmonic series, despite 
its poor convergence, the best representation for the radiation fields 
dipoles the presence cylinder whose diameter not excessively large 
compared the wavelength. With the advent automatic computing 
methods, the limitations imposed poor convergence are not stringent 
they once might have been. 

the purpose this paper present set radiation patterns 
stub antenna radial dipole located cylindrical structure. The case 
the usual isolated infinite cylinder perfect conductivity chosen first and 
compared some previously reported experimental measurements (Bain 
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1953). The results are then extended pair diametrically opposed stubs 
located the cylinder. Following this, the case half-plane sheet with 
cylindrical tip considered where the stub the cylindrical surface. 
Whereas the stub isolated cylinder considered idealized 
representation for radial antenna aircraft fuselage the stub the 
cylindrically tipped half-plane considered have some correspondence 
with antenna located near the leading edge aircraft wing vertical 
stabilizer. 

The mathematical derivation the formulae used outlined the appendix. 


DISCUSSION THE CALCULATED PATTERNS 

The explicit formulae for the radiation far fields electric dipole 
located the presence perfectly conducting circular cylinder infinite 
length were first given Carter (1943). His expressions were obtained 
application the reciprocity theorem the classical solution the diffraction 
plane wave arbitrary incidence circular cylinder. results 
have also been obtained Lucke (1951) who used different method. 

Choosing cylindrical coordinate system the circular cylinder 
occupies the space The stub antenna considered approximated 
infinitesimally short electric dipole length located (a, 
and oriented radially. The radiation field this particular antenna is, 
general, elliptically polarized. the principal equatorial plane 0), 
the radiation field plane polarized and the electric field has only compo- 


nent and proportional the factor given 
where the Hankel function the second kind integral order and 
where the argument given 


with being the free space wavelength and the diameter the cylinder. 
The function was calculated the Ferranti computer for values 
0.199, 0.800, 2.09, 5.74, and 9.67 corresponding diameter wavelength 
ratios 0.063, 0.25, 0.67, 1.83, and 3.08 respectively. ranged from +90° 
—90° intervals 5°. The choice values was dictated partly the 
sizes cylinder employed earlier experimental investigation. The 
for the five values These are definition the radiation patterns 
the stub antenna the isolated circular cylinder. Fig. the direction 
vertically upward indicated the sketch the stub antenna 
the cylinder the bottom the figure. 

observed that the pattern Fig. for 0.063 almost identical 
shape the function which the pattern isolated electric 
dipole the plane the dipole. This fact also evident examining equation 
(1) for under the condition that which case only the first term 


a. 
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(20° 


210° 30° 


Calculated 
Experimental (Ban) 


Fic. radiation pattern single stub antenna the principal plane (transverse 
the axis the cylinder) showing both calculated and measured data. (The quantity plotted 
proportional where the azimuthal angle.) 
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the summation (i.e. significant. can also seen from Fig. 
that the cylinder increases size the pattern changes its shape although 
for 0.25 the distortion not marked. For larger cylinders the field 
the rear becomes considerably less magnitude with corresponding 
increase the number lobes. interesting note that the rear lobes for 
the largest cylinder are intervals approximately 18° which corresponds 
approximately one half-wavelength arc the circumference the 
cylinder. This would expected the waves are considered travel around 
the cylinder with phase velocity near that plane wave free space and 
form interference standing wave pattern. 

also interesting note that the pattern the forward direction (up- 
ward Fig. for the larger cylinder again almost identical the cos 
pattern. This would expected physical grounds since the front face 
the cylinder would behave more like flat ground plane the radius 
curvature becomes larger compared the wavelength. 

Measured patterns for the cylinders 0.25, 0.67, 1.83, and 3.08 
that were reported Bain (1953) are shown dashed lines Fig. along 
with the computed curves. The experimental and theoretical curves are 
matched that the maximum fields are the same. The agreement quite 
reasonable. The small discrepancy can probably attributed partly the 
fact that Bain’s experiment the stub was length whereas the calcu- 
lations the stub was considered very short compared the diameter 
the cylinder. Remembering that free space, the \/4 thin monopole 
flat ground plane characterized the pattern 


sin 

whereas the infinitesimal monopole characterized cos would 
expected that the experimental patterns would slightly more directional 
towards the tangent plane 0°). This evidence the bottom two 
patterns Fig. course, would have been possible have rigorously 
considered the \/4 stub the theoretical analysis. This, however, would lead 
considerably more complicated formula for computation, and furthermore 
some particular type current distribution has assumed along the 
antenna. 

The situation where two stubs are located the cylinder and diametrically 
opposed shown Fig. this case the currents the stub are con- 
sidered the same magnitude and phase. The field the principal 
quantity specifies the pattern and shown plotted Fig. for the same range 
Fig. interesting note that the pattern for the smallest 
cylinder almost identical the function cos which characteristic 
the radiation field electric quadripole. the cylinder increases size, 
the patterns become more complicated shape, but they retain the essential 
feature that the field zero the four principal directions (i.e. 0°, 90°, 
180°, and 270°). 
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150° 120° 90° 60° 


240° 270° 


Fic. The radiation patterns the principal plane pair diametrically opposed 
stubs the cylinder fed phase. 
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\ 


Fic. radiation pattern the principal plane stub semicylindrical boss 
located flat conducting ground. (These patterns are also applicable pair diametri- 
cally opposed stubs isolated cylinder fed antiphase.) 


150° 120° 90° 6° 
180° 0° 
150° 120° 60° x 
150° 120° 90° 60° x” 
| 
3.08 
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Another very interesting case indicated Fig. where the stub situated 
semicylindrical boss located perfectly conducting flat ground plane. 
The pattern the principal plane now proportional the magnitude 
the function which shown plotted Fig. for five values 
the diameter the cylindrical boss. The formula for the radiation field 
the stub the cylindrical surface can obtained special case the 
general result derived the appendix for radial dipole located cylindri- 
cally tipped wedge where the wedge angle simply 180°. course, this 
special case also equivalent pair diametrically opposed stubs located 
the isolated cylinder and fed antiphase with equal currents. This con- 
figuration indicated the bottom Fig. and the patterns for this case 
include the dashed portions the curves. 

For small cylinders the patterns Fig. are seen again closely 
represented the function cos which characteristic short monopole 
flat ground plane. the cylinder increases size the lobe structure 
becomes more complicated, but retains the important feature that the field 
has maximum the broadside direction. 

Another case considered where the stub located cylinder which 
itself coaxial with the edge perfectly conducting half-plane. The expres- 
sion for the radiation field clearly also special case the general wedge 
solution given the appendix where the angle the wedge 360°. The 
pattern then proportional the magnitude the quantity which 
given 


with being the angular location the stub the cylinder measured from 
the half-plane and the azimuthal angle which also measured from the 
half-plane. The summation over integral values the Hankel functions 
are both integral and half-integral order. 

Computations for the function were also carried out the Ferranti 
computer for the two cases where 90° and 180°. The patterns for these 
are shown Figs. and noted that for the smallest cylinder, the 
pattern Fig. closely represented the function which has 
the angular dependence monopole the edge half-plane sheet 
(Johnson 1955; Wait 1955). This fact also evident from equation (2) when 
the first term only considered. the cylinder increases size, lobes appear 
the direction the rear. striking that the patterns the forward 
direction (to the left Fig. are characterized closely the pattern the 
stub isolated cylinder. That is, the effect the half-plane only 
noticeable for directions the rear (to the right Fig. 4). 

Fig. the stub located 90° from the half-plane and the patterns are 
changed radically shape. this case there discontinuity the field 
along the half-plane sheet. For small cylinders, the pattern beginning 
approach cos(¢/2) function, but the effect the asymmetry the location 
the stub with respect the half-plane still quite noticeable. For larger 


(2) 
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240° 270° sor 210° 


Fic. The radiation pattern the principal plane stub the cylindrical surface 
tipped conducting half-plane. (The stub located 180° from the half-plane.) 
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210° 240° 270° 300° 330° 


Vhe radiation pattern the principal plane stub the evlindrical surface 
conducting stub located from the half-plane.) 
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cylinders, multiple lobes appear only the upper right-hand quadrant. 
fact, evident that the pattern this quadrant very similar the 
patterns shown Fig. for the stub the semicylindrical boss, flat 
ground plane. the other hand, the pattern the other three quadrants 
for the larger cylinders quite similar the patterns Fig. for the stub 
the isolated cylinder. 


CONCLUDING REMARKS 


The patterns Figs. should prove useful making rough estimates 
the patterns for obstacles more complicated shape. The absolute magni- 
tude the radiated fields has not been considered, since this usually second- 
ary interest. The numerical values the functions 
and are file the Computation Centre the University Toronto. 
would interest extend these calculations other than the principal 
plane where the field is, general, elliptically polarized. this case would 
necessary specify not only the magnitude the two cross-polarized com- 
ponents the field, but also the relative phase between them. 
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MATHEMATICAL APPENDIX 

The general problem calculate the fields (see Fig. Al) due 
element current with respect conducting wedge exterior angle 
with cylindrical tip radius Choosing cylindrical coordinate system, 


\ 
Nhe 


Fic. Al. Wedge with cylindrical tip. 


the surfaces the wedge are defined and and the tip repre- 
and the tipped wedge reduces semi-infinite sheet. instructive 
let the magnitude and remain arbitrary for the general analysis. 
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Following previous analyses (Wait and Kahana 1954; Wait 1955) the axial 
components the electric field and the magnetic field, the region bounded 
the wedge, can expressed by* 


V € 
and 
Vv 


where and are vectors extending from the origin and respectively 
and where and p(r’) are the electric current density vector and the 
electric charge the point P’. The gradient operator involves differentia- 
tion with respect the primed coordinates. The Green’s functions r’) 
and r’) satisfy the inhomogeneous Helmholtz equation 
where 7’) the three-dimensional Dirac function. The functions 
are subject the following boundary conditions: 


the above, the time factor and the medium contained the 
volume homogeneous with electrical constants and and wave number 
space wavelength). 
The appropriate forms the Green’s functions for wedge regions were also 
given previously (Wait and Kahana 1954). the far-field region, such that 


and 


_ J,’ (ua) ] 
function the first type, and the Hankel function the second 
kind. 


*N.B. should pointed out that the earlier paper (Wait and Kahana 1954), the 
currents were magnetic type. 
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When these expressions for and are substituted equations (3) and 
(4), the complete far-field solution obtained for any distribution current 
within the wedge region. The situation, illustrated Fig. A2, now specialized 
source consisting thin linear antenna carrying current con- 


Fic. A2. Linear antenna radial tip wedge. 


angle The volume integrals over then become line integrals over 
the length the antenna. then follows without difficulty from equations 
(3) and (4) that 


and 


where the equation continuity relating charge and electric current has been 
utilized. The field components the far zone, terms spherical coordinates 
(R, are given 


The final results for the radiation fields are then 


and 
. —ikR Po on 
1 m= 
where 
and 
(14) G,(u) Ee )-H, (up ) (ua) up’ 


P, 
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The above formulae are quite suitable for computation p2, and are 
less than wavelength so, since then the infinite series converges very 
rapidly. fact, p2, and are small compared the wavelength, equations 
(11) and (12) are given adequately considering only the first term; then 


and 
where 
e 


independent 

The general expressions for and G,(u) given above equations (13) 
and (14) can considerably the length the antenna can 
considered short compared the radius the cylindrical surface. this 
case p’, where occurs equations (11) and (12), can replaced 
that 

where Wronskian identity has been utilized. The corresponding change 
the principal plane 90°), the fields are now given 


(18) 


(19) 


where ka. Note that for the case semicylindrical boss flat ground 
and m/2. 
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RANGE-ENERGY RELATIONS FOR LIGHT NUCLEI 
AIR AND NUCLEAR 


ABSTRACT 


empirical method has been developed for estimating the effective charge 
nuclear emulsion. The results are used conjunction with power-law function 
for the variation stopping-power with velocity derive range-energy relations 
for these ions. The curves are compared with all available experimental data, 
and agreement within the limits experimental error obtained all cases 
except those fluorine and neon ions air. 


INTRODUCTION 


The practical problem estimating the energy nuclear particle from 
its observed range known medium important that many experimental 
and theoretical studies have been made with the object establishing standard 
range-energy relations. Most these investigations referred protons and 
alpha-particles the common gases, metal foils, and nuclear emulsions 
the concentrated type. sufficiently high energies the rate energy loss 
any charged particle can described fairly accurately theoretical 
formulae employing unique value the charge. low energies, however, 
the effective charge decreased the occasional capture orbital electrons 
and adequate theoretical treatment has been given for this region. This 
paper describes empirical method for estimating the effective charge 
light nucleus function its velocity. The results are used calculate 
range-energy relations for ions nuclear charge between (lithium) 
and (neon) air and C.2 emulsion. Capture electrons into 
the and shells only considered; for ions with higher than the 
present method less valuable and the method Knipp and Teller (1941), 
based the atomic charge distribution, should preferable. 
This latter method has previously been used Wilkins (1951) and Lonchamp 
(1953) calculate range-energy relations emulsion and will interest 
compare their results with those obtained here entirely different 
method. 

believed that charged particle traversing matter loses energy chiefly 
electronic collisions and some extent nuclear collisions, although the 
nuclear effect importance only for heavy ions such fission fragments. 
Allison and Warshaw (1953) have calculated that even for neon ions 0.4 
Mev. argon the loss nuclear collisions only 17% the total energy 
loss; accordingly these effects are not considered separately this paper, 
view the considerable errors available experimental data the subject. 
The various theoretical formulae for the rate energy loss charged particles 
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electronic collisions have been reviewed Uehling (1954) and the condi- 
tions validity for each type formula have been studied Bohr (1948). 
general one may write for given medium: 


where the energy loss per unit distance traversed, 

the effective charge the particle, 
and F(V) some function the velocity depending the mean ionization 
potential the medium (J). 


For comparison with experimental data assumed that the net charge 
ion equilibrium given velocity with surrounding electrons, but 
this does not necessarily correspond the actual conditions along the track 
particle being slowed down dense medium. 

Protons have unit charge over wide range velocities and experiments 
with these particles determine F(V) with some accuracy down the region 
where the same order magnitude the classical orbital velocity 
the electron the hydrogen atom (Vo). According Bohr (1940) capture 
and loss electrons should equally probable when close that is, 
should approximately 0.5 when the proton energy 0.025 Mev. Experi- 
mental the rates energy loss and the effective charges different 
ions have become widely available the review article Allison and Warshaw 
(1953). Their figures are used extensively ensuing sections this paper and 
theoretical predictions concerning effective charges are compared with the 
estimates Gluckstern (1955) who used the method Bell (1953). 

With reliable information about the behavior function velocity 
possible integrate the reciprocal energy loss dx/dE over the appro- 
priate energy intervals obtain relation. The results such 
calculations for air are compared with the available experimental data 
Section III, and those for Ilford C.2 emulsion Section IV. accompanying 
paper the identification (y, stars emulsions illustrates the 
practical applications curves for lithium, beryllium, and boron 
ions. 


THE BEHAVIOR CHARGED PARTICLES AIR 


Because the difficulties attendant upon the full treatment energy losses 
slow charged particles has been found necessary adopt semiempirical 
method based equation (1). For protons air 15°C. and standard 
pressure the function approximates power-law the form: 


(2) F(V) 


dE/dx Mev./cm., and 1.34. The choice not entirely arbitrary 
because its value may checked against the formula Bethe (1930) which 
believed hold for fast protons well above the charge-exchange region. 
Direct analysis the energy-loss data Allison and Warshaw (1953) and the 
relation Bethe and Ashkin (1953) gave 1.34, while the 
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Bethe formula with ev. gave 1.36 for protons between 0.5 Mev. 
and 2.0 Mev. The specific stopping cross-section used Knipp and Teller 
(1941) approximates power-law with 1.36, although very low 
velocities the onset nuclear collisions should lead higher value 

The really uncertain point the use equation (2) lies the extrapolation 
the power-law relation low velocities where charge-exchange important 
and the Bethe formula longer valid. Here essential have check 
the method and this provided the experimental measurements the 
effective charge low velocities. Having chosen 1.34 can use the 
function extrapolated downwards together with the known energy-loss 
data find the value any velocity. Fig. the calculated values 
are plotted against and the smooth curve may compared 
with the experimental points quoted Allison and Warshaw (1953) and 
Stier, Barnett, and Evans (1954). clear that the main trend the curve 


Fig. 1A. PROTONS IN AIR 
calculated 
curve 
05 Experimental points 


70° e Stier, Barnett and Evans 
Allison and Warshaw 
Fig. 1B. HELIUM AND NEON IONS AIR 
= 
calculated 
helium 
curve Experimental points 
Stier, Barnett and Evans 
Helium Allison and Warshaw 
os Knipp and Teller 
Stier, Barnett and Evans 
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Particle velocity Vv 


Fic. Variation effective charge with velocity calculated from energy-loss data and 
compared with experimental points. Protons air. Helium and neon ions air. 
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curve 


206 CANADIAN JOURNAL OF PHYSICS. VOL. 34, 1956 


correct down values below unity and the existing discrepancies 
can accounted for postulating contributions from nuclear collisions 
since higher value would improve the agreement. The smooth curve has 
been extrapolated zero such fashion that trial reintegration the 
proton curve gave results consistent with observed values. 

deduce from the foregoing analysis that the effective charge ratio 
related functionally the velocity ratio V/Vo way agreeing 
essentially with the use equation (2) although detailed agreement not 
attained. When extend the idea using such curve the case helium 
ions, however, have two orbital electrons take into account, and charge- 
exchange important even particle energies around Mev. Throughout 
this work the particle velocity related which the velocity the 
first electron the appropriate shell, because practice approach 
the electron-capture problem from the high-energy end where this signi- 
When the analysis repeated for helium ions, with 1.34 before 
equal twice the value used for protons, the calculated Z,/Z curve lies 
well above the proton curve because the second electron has lower velocity 
than the assumed Fig. the smooth curve may compared with 
experimental values quoted Knipp and Teller (1941) and Allison and 
Warshaw (1953), also with the recent data Stier al. (1954). general the 
calculated curve behaves correctly but everywhere higher than the experi- 
mental points; this effect may again attributed nuclear collisions which 
are more important here than the proton case. higher value would 
give better agreement but was decided retain 1.34 and maintain 
consistent scheme calculation. 

The chief difficulty using the curves Fig. for deriving 
relations lies the unknown behavior ions very low velocities. The 
theoretical work Fermi and Teller (1947) shows that the power-law equation 

cannot valid all very close the origin Fig. although, will 
seen later, necessary cut off the divergent integral low 
energies that all errors this region are combined into one additive error. 
The comparison with the experimental results rendered uncertain the 
fact that the data Stier a/. (1954) differ markedly from those Snitzer 
(1953) the low-energy region. Moreover, more and more electrons are 
captured into shell, their average velocity reduced and the inner shells 
are also affected. The chosen value therefore not reliable guide, 
may seen from the neon data included Fig. 1B. Here the observed values 
Z,/Z for neon ions are plotted against with corresponding the 
ionization potential the first electron the shell (238 ev.). These points 
lie well the left the helium points although both sets refer full electron 
shells. order able extend the scheme calculation ions with 
L-shell capture, the energy losses neon ions quoted Allison and Warshaw 
(1953) were used calculate curve for near the origin; the line obtained 
agrees quite well with the experimental points Stier (1954) and 
joins the helium line satisfactorily. The final line lies between the sets 
points for helium and neon ions and subsequent calculations this has been 
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taken typical curve for full shell, conjunction with the 
equation (2). 

order check the method calculating very low energies 
the energy loss nitrogen ions air was estimated and the results are shown 
Table together with the experimental figures quoted Allison and War- 


TABLE 
RATE ENERGY LOSS NITROGEN IONS AIR 
calculated from equations and 


Mev. energy loss energy loss 
0.20 1.66 0.7 2.3 2.83 
0.30 0.9 3.0 3.47 
0.40 2.34 1.0 3.3 4.00 


shaw (1953). this work the possible ionization the shell was neglected 
and curve intermediate between the curves Fig. and was adopted 
for the ionization since only five the eight electrons are available 
nitrogen. view the extreme difficulty estimating the effective charge 
this region the agreement between theoretical and experimental values 
must regarded satisfactory. small contribution for A-shell 
ionization would make the agreement very close. 
RANGE-ENERGY CALCULATIONS FOR AIR 

calculating the range-energy relation for any particular ion air, the 
velocity the first orbital electron the shell was found the 
ordinary Bohr formula; the value for the first electron the shell 
was estimated from the ionization potential figures Landolt and 
Tables (1950). The velocity the particle was expressed fraction 
these two values with extra allowance the case for small 
increase when the became appreciably ionized. This introduced 
correction the final results but was importance only the curves for 
The effective charge the shell, was found from the 
smooth curve Fig. and the value for the shell, was estimated 
interpolation between the curves Fig. and according the degree 
occupation the shell the neutral ion. Finally was calculated for 
the whole ion and was tabulated function with the aid 
equation (1) and (2). 

the final integration over the entire energy range from zero 
given value was necessary deal with the divergence the 
expression dx/dE approaches zero. standard procedure for cutting off 
the integral linearly was worked out reintegration the curve for helium 
ions and this procedure gave results very nearly independent the cutting-off 
point provided that this was chosen the region where dx/dE reaches twice 
its minimum value. The error introduced here additively constant for any 
one value and believed less than the corresponding experimental 
errors most cases. 


| 
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The range-energy relations for ions between and were 
finally expressed form independent the isotopic mass plotting 
compared with the results Knipp and Teller (1941) the modified theory 
Brunings, Knipp, and Teller (1941) and found that they agree fairly 


RANGE ENERGY RELATIONS FOR IONS Z=3 AIR 
15°C, 
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Fic. relations for ions air compared with experimental points based 
cloud-chamber measurements. 


closely with the 1.1 curves the earlier paper. Experimental data for 
comparison with the theoretical curves are not plentiful and unfortunately 
all the published figures are based cloud-chamber measurements which 
stopping-power suitable for fast particles used convert observed ranges 
air ranges’’. Fig. shows series experimental points for 
0.1 Mev. obtained for carbon ions Feather (1933) and Wrenshall 
(1940), for nitrogen and oxygen ions Blackett and Lees (1932), for fluorine 
Feather (1933), and for neon Eaton (1935) and McCarthy (1938). 
The agreement good for the data referring ions with and but 
marked discrepancies occur for the fluorine and neon results. The explanation 
may lie either the use the incorrect stopping-power conversion the 
ranges air, the possibility that the filling-up the last electrons the 
shell proceeds more slowly than indicated Fig. 1B, although the neon 
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results low energies were incorporated the curve effort make 
reliable possible for L-shell calculations. Typical values for 
below 1.0 Mev. are shown Table 


TABLE 


CALCULATED RELATIONS FOR IONS AIR 
15°C. and standard pressure 


E/M, cm. for ions with 
Mev 
0.1 0.28 0.39 0.55 0.77 0.95 1.16 1.9 
0.15 0.40 0.53 0.75 1.00 1.25 1.52 1.9 2.4 
0.25 0.63 0.82 1.10 1.45 1.80 2.15 
0.5 1.28 1.6 2.0 3.0 4.3 
0.75 2.0 3.0 3.6 4.2 4.9 6.7 
1.0 2.9 4.1 4.8 5.5 6.3 8.3 


IV. RANGE-ENERGY RELATIONS NUCLEAR EMULSION 


the method the foregoing sections used for 
calculations emulsion, must know the appropriate function F(V) and 
also have some information about the charge-exchange properties ions 
the medium. Although extensive measurements have been made the ranges 
fast protons emulsions and theoretical calculations are generally 
good agreement with the results, low velocities the behavior protons 
such inhomogeneous medium cannot safely predicted. further 
culty the use existing curves for helium ions that these 
have usually been constructed from series measurements scattered 
energy intervals with the aid the assumption that the integral stopping- 
power varies smoothly with energy. Below Mev. such curves are limited 
value for analysis. additional difficulty arises from the fact that deter- 
minations effective charge have been carried out emulsion and the most 
numerous data available refer metal foils, which the process electron- 
capture probably differs greatly from that gelatin and silver bromide. 

The proton relations Wilkins (1951) and Gibson, Prowse, 
and Rotblat (1954) were used find F(V) the form power law and this 
agreed with equation (2) when the constants were chosen follows: 1.15, 
and 1.04 units when was units cm./sec. and dE/dx was 
Mev./micron. The formula Bethe (1930) with 400 ev. gives values 
between 1.0 and 1.1 the appropriate energy region. Various measurements 
the ranges alpha-particles having energies between 2.1 Mev. (samarium) 
and 8.78 Mev. (thorium C’) were found fit empirical equation the 
form: 


and assume that alpha-particles are completely stripped electrons above 
Mev. this leads the relation: 


—dE/dx 


: 
| 
| 
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This value rather higher than that for protons, but, the calculations 
for air, chose the proton figure finding the curves for protons and 
helium ions emulsions. Fig. the results for 1.15 may compared 
with the experimental results for metal foils quoted Allison and Warshaw 
(1953) and Dissanaike (1953). Near the origin the figures for neon ions 


PROTONS, HELIUM AND NEON IONS EMULSIONS 


calculated 


calculated proton curve 


helium 
curve 


Experimental points (metal foils) 


Protons @ Allison and Warshaw 


Allison and Warshaw 
Helium 
Dissanaike 


effective charge ratio 


Particle velocity 
Orbital velocity 


Fic. effective charge with velocity for protons and helium ions emulsion 
compared with experimental points for the same particles metal foils. 


argon gas have been included the grounds that these ions seem carry 
the same charge the same velocity several different media. This point 
has been studied specifically Hubbard and Lauer (1955) who found 
essentially the same stripping behavior oxygen and neon ions several 
gases. the same time one would expect condensed media behave rather 
differently, according Bohr and Lindhard (1954), and metals may differ 
from non-conductors considerably. Fig. there are obvious discrepancies 
between the curves calculated for emulsion and the experimental points 
obtained with metal foils. 

order investigate the reliability Fig. for estimating effective 
charges the shell, calculations were first made for nitrogen ions various 
velocities and Table III the results are compared with the most extensive 


TABLE 


CALCULATED EFFECTIVE CHARGE NITROGEN IONS 
Compared with the results Reynolds al. (1954) for nickel 


Energy Mev. 1.2 2.4 


Velocity 
Calculated 3.3 4.4 4.8 5.3 5.7 6.1 6.3 6.4 
Observed (2.0) (8.2) (4.2) 4.9 5.3 5.5 5.7 5.9 6.2 


10 
O 
a 
e 
a 
Neon ¢ + 
NIN curve , 
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set figures present available, namely those Reynolds, Scott, and 
Zucker (1954) nitrogen ions nickel. The low-energy values are based 
rather uncertain extrapolation the range—velocity curve and these are 
enclosed brackets. The results show systematic error giving too high 
calculated values both high and low velocities although the agreement 
good intermediate velocities. 

attempt find out the discrepancies are due the use data 
obtained with metal foils, separate set calculations was made for nitrogen, 
oxygen, fluorine, and neon ions and these are compared Table with 


TABLE 
EFFECTIVE CHARGE OF VARIOUS IONS IN GASES AND ORGANIC FILMS 
4.7 4.5 4.0 4.68 (a) 
5.4 4.9 §.33 (a) 
12.7 5.6 5.5 5.35 
§.7 5.6 5.40 
5.9 5.8 5.6 5.49 (a), 
6.2 6.1 6.09 (a) 
6.4 6.3 6.19 (a) 
6.1 4.7 4.8 4.4 4.8, 4.86 (c), 
8.7 5.4 §.1 5.2 (c) 
17.1 6.5 6.4 6.22 
20.4 6.9 6.63 
§.1 5.4 4.7 4.9 (c) 
10.3 5.9 6.0 6.0 (c) 


(b) Stephens and Walker 1955. 
(c) Hubbard and Lauer 1955. 


experimental data referring ions emerging from organic films and travelling 
through gases. These data were obtained Stephens and Walker (1955), 
Reynolds, Wyly, and Zucker (1955), and Hubbard and Lauer (1955). 
the table Z,(air) gives the results calculations made with the curves 
Fig. and Z,(em) those from Fig. while values are the theoretical 
predictions Gluckstern (1955) for ions argon. will seen that the 
present theory gives rather too high values for and that the results 
Gluckstern are general too low, but nevertheless the agreement satisfac- 
tory. would appear that the curves Fig. are reliable for estimating 
L-shell ionization with accuracy the order 10% despite the discre- 
pancies between the curves and the data for metal foils. 

Regarding the foregoing results sufficiently accurate allow the calcula- 
tion relations for emulsion, have repeated the work 
Section III for this medium. Trial reintegrations the proton and helium 
curves gave figures which were correct within over the energy range 
between E/M 0.1 Mev. and E/M 1.0 Mev. The results for the heavier 
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~ 


ions are shown Fig. Particular attention should paid the case 
because the work Reynolds and Zucker (1954) has yielded 
curve for nitrogen ions Mev. The two curves shown for 
were obtained: with small correction for the change Vo” the K-shell 
ionization becomes more pronounced, and without this correction, 


RANGE ENERGY RELATIONS FOR IONS ILFORD 


EMULSION 


Experimental points 


in microns 


Ol o2 03 04 05 06 0-8 


in Mev. 


Fic. Range-energy relations for ions Ilford C.2 emulsion. The experimental points 
Reynolds and Zucker (1954) for nitrogen ions may compared with curves for 
(calculated with correction for K-shell ionization) and (without the correction). 


being fixed this case the ionization potential the first electron 
(98 ev.). The experimental points lie close these curves and the discrepancies 
are within the limits experimental error values E/M near 2.0 Mev. 

Another accurate comparison with experimental data afforded the 
observed ranges lithium nuclei produced the reaction: 


Fig. the accompanying paper the results are compared with other 
experimental determinations and with the theoretical line. low energies 
there small discrepancy between this theory and experiment and the 
curve Wilkins (1951) seems slightly more accurate but the differences 
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are the order 0.5 micron only observed range. The present results 
differ markedly from those Locher and Stoll (1953) who seem have used 
constant ionic charge calculating their theoretical curve. 

comparison the present theory with that Wilkins (1951) shown 
for the case carbon ions Fig. this paper, where the experimental 
points are chiefly due Miller (1952) with two extra points fixed large- 


CARBON IONS IN EMULSIONS 


30 
Experimenta! points 
aa From inelastic proton scattering 
Miller 
3 20 
N 
Theoretical lines 
present theory 
‘9 ----— Miller 


—-— Wilkins 


=z in Mev 


Fic. relations for carbon ions emulsion. Theoretical lines from Miller 
(1952), Wilkins (1951), and this paper may compared with experimental points due 
Miller and two cases inelastic proton scattering emulsion. 


angle scattering protons emulsion observed the author. The calculated 
line Miller (1952) also shown and the work Lonchamp (1953) gives 
line close that Wilkins. the whole the present theory fits the results 
better than the others although gives values too small low energies. 
likely that the method Knipp and Teller (1941) gives consistently low 
values for the range the energy region Mev. (J. Wilkins, 
private communication). 

further check the theoretical relations high energies provided 
the observations Barkas (1953) the range corrections applied 
values regions where electron-stripping complete. For these 
conditions one can write: 


where constant depending only the value Barkas quotes values 
for and and these are compared with the predictions this 
theory and that Lonchamp (1953) Table The present theory gives 
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TABLE 
RANGE CORRECTION Cz FOR IONS AT HIGH ENERGY 
Z?R/M, values microns 


(this paper) (Lonchamp) Observed value 
3.4 2.7 

8.4 6.2 

15.4 10.5 13.5 

26.5 15.9 23.9 


rather high values but they agree with experiment within micron 
the observed range, and this considered reasonable under existing 
experimental conditions. 

Typical values different energies are quoted Table VI, 
and these are believed substantially accurate for ions 


TABLE 
CALCULATED RANGE-ENERGY RELATIONS FOR IONS EMULSION 
0.1 3.0 4.4 5.9 8.9 10.7 
0.15 2.9 4.0 5.0 6.0 9.6 11.7 13.8 
0.25 4.4 5.8 7.2 8.7 11.0 13.4 16.0 18.8 
0.5 8.7 10.6 12.8 18.3 21.8 25.5 29.5 
0.75 13.3 18.4 21.8 25.9 29.6 33.7 38.5 
1.0 18.6 21.3 24.5 28.5 33.2 37.6 42.3 47.5 


Unfortunately have reliable data for fluorine and neon ions emulsion 
light can thrown the discrepancies between theory and experiment 
noted the case air. 


CONCLUSION 


The theory used calculating the rates energy loss various ions air 
and nuclear emulsion assumes power-law velocity function and definite 
relation between the effective charge given electron shell and the ion 
velocity expressed fraction the classical velocity the first orbital 
electron. Both these assumptions are subject considerable uncertainty 
but series comparisons with experimental data shows that the method 
least accurate other methods when applied light nuclei below 
10. Various calculations using different methods now agree fairly closely 
about the form relations for ions but there are 
difficulties accounting for the ranges fluorine and neon ions air, and 
set experiments determine these would particular value. 


ACKNOWLEDGMENTS 


wish thank Dr. Sargent and Mr. Allen for valuable discussions 
and Dr. Wilkins for private correspondence about the experimental data. 
indebted the National Research Council for financial assistance. 


to 


LIVESEY: RANGE-ENERGY RELATIONS 


REFERENCES 


W.H. 1953. Phys. Rev. 

BELL, G.I. 1953. Phys. Rev. 

H.A. 1930. Ann. Physik, 

and AsHKIN, 1953. Experimental nuclear physics. Vol. Edited 
Segré. John Wiley Sons, Inc., New York. 

M.S.and LEEs, D.S. 1932. Proc. Roy. Soc. (London), 134: 658. 

1948. Danske Math.-fys. Medd. 18, No. 

1954. Danske Math.-fys. Medd. 28, No. 

M., K.,and TELLER, 1941. Phys. Rev. 60: 657. 

Eaton, 1935. Phys. Rev. 

FEATHER, 1933. Proc. Roy. Soc. (London), 141: 194. 

E.and TELLER, 1947. Phys. Rev. 

GLUCKSTERN, R.L. 1955. Phys. Rev. 

Knipp, K.and TELLER, 1941. Phys. Rev. 59: 659. 

1950. Zahlenwerte und Funktionen, Auflage, Band—Atom und 
Molekular Physik, Teil—Atome und Ionen. Verlag von Julius Springer, Berlin. 

1953. Phys. Rev. 90: 164. 

1938. Phys. Rev. 53:30. 

MILLER, 1952. Univ. Calif. (Berkeley) Rad. Lab. Rept. UCRL 1902. 

Phys. Rev. 89: 1237. 

STEPHENS, G.and WALKER, 1955. Proc. Roy. Soc. (London), 229: 376. 

UEHLING, 1954. Ann. Rev. Nuclear Science, No. 

WILKINS, J.J. 1951. A.E.R.E. Rept. G/R 664. Harwell. 

WRENSHALL, 1940. Phys. Rev. 57: 1095. 


REACTIONS PRODUCED 
PHOTOGRAPHIC EMULSIONS! 


ABSTRACT 
Ilford C.2 plates were exposed bremsstrahlung maximum energy Mev. 

from the electron synchrotron Queen’s University, and stars produced 
the emulsion have been attributed the reaction: 

Mev. 
The excitation function the reaction shows sharp drop near Mev. and 
the integrated cross-section this limit has been estimated. addition 
stars have been identified being due the reaction: 

Mev. 
both reactions the most probable mode disintegration the initial 
emission alpha-particle. The range measurements the carbon stars have 
been used derive accurate range-energy relation for lithium ions nuclear 
emulsion. 


INTRODUCTION 


Among the many types charged-particle stars produced exposing 
photographic emulsions high-energy X-rays one shown schematically 
Fig. consisting proton track, alpha-particle track, and heavier 


TYPICAL STAR 


recoiling light 
nucleus eg LI 


y momenta 
proton 


x momenta 
perpendicular 


‘ 
recoil 


Fic. Diagram show the mode occurrence (y, p+a) stars and the method 
momentum balance used their identification. 


recoil. The principle identifying these stars with the aid the 
momentum relations for the recoil nuclei has been described Livesey and 
Smith (1952) and accompanying paper the present author (1956) out- 
lines empirical method for calculating the relations for light 
nuclei nuclear emulsion. The derived curves are here 
employed identify stars attributed the following reactions: 
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Silver and bromine also yield (y, stars but the recoils are always too 
short measured accurately. The other isotopes carbon, nitrogen, and 
oxygen have also been ignored the present work. Little information 
yet available about the three reactions quoted, although Brinkworth, Goward, 
and Wilkins (1952) have discussed the difficulties identifying the stars 
low energies, and Goward and Wilkins (1955) gave the inte- 
grated cross-section for the carbon reaction. The present experiment employed 
X-rays energies Mev. and yielded new results chiefly about the 
carbon and nitrogen reactions. incidental result the star analysis, the 
range-energy relation for lithium ions C.2 emulsion has been de- 
termined for energies below Mev., and the curve agrees well with the theo- 
retical relation found the accompanying paper (Livesey 1956). 


EXPERIMENTAL DETAILS 


Photographic plates measuring in. in. and coated with 200-micron 
layers Ilford C.2 emulsion were exposed bremsstrahlung maximum 
energy Mev. produced the electron synchrotron Queen’s University. 
collimated beam was used and each plate received estimated dose 
roentgens. Most the plates were placed right angles the beam 
with the emulsion side facing the synchrotron target, but number stars 
were studied two plates exposed parallel the beam. these plates the 
principle momentum balance used star identification had modified 
slightly, but the results were essentially similar those the other plates, 
and the main difference lay the angular distributions individual particles 
belonging accepted stars. was found that preference was given stars 
with tracks nearly parallel the plane the plate, whatever the direction 
the incident photons may have been. therefore difficult obtain genuine 
angular distributions unless series plates exposed several angles the 
beam are used. this experiment the number stars identified was too small 
provide reliable information this kind. 

After exposure the plates were processed method described Dawson 
and Livesey (unpublished) which gives good discrimination between the tracks 
different particles without appreciable loss fast proton tracks. The plates 
were scanned two observers using microscope magnification and 
all suitable stars were measured accurately high power. Special attention 
was paid the problem measuring short recoil tracks and difficult cases 
high-power eyepiece was used. The errors range measurements are esti- 
mated +0.5 micron except stars where the recoils dipped steeply the 
emulsion, when identification almost impossible. recoil range less than 
3.0 microns was considered the final analysis, and stars which had track 
escaping from the emulsion boundaries were not measured. 

practice the observer first identified star visually belonging the 
group and the balance momentum between the proton and alpha- 
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particle the horizontal plane was checked, under the assumption that 
neutrons were emitted the reaction. The star was accepted the end 
the analysis the net error momentum balance, compounded from the 
the and directions, did not exceed photon units, i.e. one 
unit momentum, but this preliminary stage the star was 
rejected unless the (in Fig. balanced within +20 photon 
units. The chief ambiguity lay the possibility (y, reactions occur- 
ring, and particular care was taken eliminate stars where the identification 
the light particle deuteron would give better momentum balance. 
The momentum were used find the momentum the 
recoil nucleus, and the obtained was plotted against the measured range 
graph similar Fig. Here the experimental points may compared 
with theoretical curves for nuclei and and nearly all cases 
agreement with one these curves found allow error micron 
the ranges. The energy the photon responsible for particular star was 
calculated the assumption that all three particles were created 
ground states, and this figure was used the final check the momentum 


Comparison of experimental points with theoretical curves for it 
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Fic. Typical range-momentum plot used identify the recoil nucleus. Experimental 
points are compared with theoretical curves obtained from the results the accompanying 
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balance for the third (Z) direction. previously stated, stars were rejected 
their net error momentum balance exceeded photon units. 

The points scattered around the curve Fig. have been used de- 
termine the range-energy curve more accurately than the theory allows, and 
Fig. the solid experimental line may compared with the dotted line 


20 


Experimental data 
Full line from this experiment 
© Barkas 4 Neuendorffer et al. 
Faraggi Cuer and Lonchamp 


im microns 


M 


A eoretical 
line 


in Mev. 


Fic. relation for lithium ions C.2 emulsion expressed 
(ordinates) plotted against (abscissae). 


given the theory, also with other experimental points published Far- 
aggi (1949), Barkas (1953), Ciier and Lonchamp (1951), and 
Inglis, and Hanna (1951). The various values are close agreement and the 
relation believed accurate within micron over the range energies 
investigated. 

Approximately 300 stars were measured series plates exposed the 
Mev. X-rays, and, these, stars were accepted due the carbon 
reaction, while were accepted nitrogen stars. detailed analysis 
oxygen stars was attempted with the small amount information obtained 


about them. 
III. RESULTS AND DISCUSSION 


The stars accepted due the carbon reaction were first used plot 
yield curve, i.e. the number stars per Mev. energy interval was plotted 
against the photon energy. order that this converted true cross- 
section curve, necessary apply corrections for (a) the loss low-energy 
stars owing the rejection very short recoil tracks, (b) the escape tracks 
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from the emulsion boundaries, and (c) the shape the spectrum brems- 
strahlung emitted from the synchrotron target. Fig. the numbers are 
approximately normalized 32.5 Mev. and the estimated cross-section scale 
included the right-hand ordinate. The dotted portion the histogram 


Group 


B—, 


Corbon 
stars 


35 40 


Corrected number of stors per Mev 


Nitrogen stars 


Approximote cross - section in milliborns 


25 30 35 40 45 50 
Photon energy in Mev 


Fic. Variation cross-section with photon energy for the (y, reactions carbon 
and nitrogen, based stars and stars respectively. 


shows the low-energy stars with small correction for effect (a). The same 
diagram shows the rough cross-section histogram for nitrogen stars and 
will seen that the yield this reaction low above Mev. whereas the 
carbon reaction has appreciable cross-section Mev. 

The integrated cross-sections calculated from these data are follows: 


Energy range Integrated cross-section 
Type star Mev. Mev-millibarns 


The first estimate may compared with the results Goward and Wilkins 
(1955) who quote cross-section 0.5 Mev-mb. for the energy range between 
and Mev. separate estimate the total integrated cross-section for 
the carbon reaction the highest energy recorded (48 Mev.) was obtained 
comparing the yield these stars with that stars the same 
area the photographic plate. The result was 1.6 Mev-mb., agreement 
with the above figures. 

The distributions the individual particle energies were also studied, and 
for this purpose the carbon stars were divided into energy groups shown 
Fig. although this division does not necessarily correspond with any genuine 
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change reaction mechanism. The energies obtained varied considerably 
within the groups but the majority values lay inside the following limits: 


Group Proton ion energy 
energy 
3.0+1.2 1.8+0.5 
5.0+2.0 4.2+1.5 


noteworthy that the alpha-particle usually gains more energy than either 
the other two particles, and clear from these results that the reaction 
cannot proceed direct three-body break-up except very small number 
cases. seems likely from the above evidence that the alpha-particle 
often emitted first, leaving the residual Be® nucleus excited state which 
can lead emission proton. the same time there are two alternative 
mechanisms for the carbon reaction, namely via the unstable nucleus 
via various excited states the nucleus 

order investigate these possibilities more directly the calculated exci- 
tation energies for intermediate and nuclei were plotted, 
Fig. and the histograms were compared with the positions levels known 
emit the appropriate kind particle. For example, there are levels 
known decay proton emission 0.4, 0.9, 1.9, 4.2, and possibly 5.3 Mev. 
above the threshold for this type decay 17.2 Mev. (see Ajzenberg and 
Lauritsen 1955). will seen from Fig. that there are considerable 
culties preventing unambiguous determination the mode decay 
these stars; many the stars can accounted for several ways, and 
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Fic. Histograms the calculated excitation energies intermediate and 
nuclei the carbon reaction, compared with known particle-emitting levels these nuclei. 
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significant peak one distribution likely lead non-genuine peaks 
the other two. the case however, there are two peaks the distri- 
bution, and statistical grounds one would expect least one these 
significant. Further analysis shows that the stars group the yield curve 
are responsible for the peak near Mev. the histogram, while 
the stars group contribute both the lower peak the diagram and 
the main peak the diagram. deduce that many, but not all, 
the carbon stars disintegrate via levels Mev. and Mev. 
above the threshold 17.2 Mev. 

The nitrogen stars were analyzed the same way but high statistical 
accuracy was not available and the only clear conclusion was that known levels 
could account for nearly all the observed stars. 

IV. CONCLUSIONS 

The experimental results previously quoted may used draw the 
lowing conclusions about this type reaction: 

(i) the range-momentum plot enables one identify the carbon and 
nitrogen reactions with considerable accuracy energies above Mev. 
and Mev. respectively, 

incidental result the range-energy relation for lithium ions C.2 
emulsion has been determined, 

(iii) the excitation function for the carbon reaction shows some structure 
high energies, notably sharp drop Mev., while the nitrogen 
reaction occurs chiefly limited energy region close 27.5 Mev., 

(iv) the most probable mode disintegration the carbon stars via 
various proton-emitting levels the intermediate nucleus, and the 
nitrogen stars are best explained being formed via the intermediate 
nucleus 

ACKNOWLEDGMENTS 

must thank the National Research Council for the financial assistance 
which made this work possible and Dr. Sargent for the opportunity 
work Queen’s University. very grateful Mr. Dawson, Dr. 
Bigham, and Mr. Rock for making some their unpublished 
results available, and Mrs. Edith MacPhail for measuring several stars. 
Mr. Janzen kindly operated the synchrotron for the plate exposures and 
Mr. Dawson contributed much valuable advice. This work has been part 
the synchrotron program receiving financial support from the Atomic Energy 
Control Board. 


REFERENCES 


AJZENBERG, and 1955. Revs. Mod. Phys. 

W.H. 1953. Phys. Rev. 

Harweil. 

Dawson, and Livesey, Unpublished. 

1949. Compt. rend. 229: 1223. 

D.L. 1956. Can. Phys. 34: 203. 

and 1952. Proc. Phys. Soc. (London), 65: 297. 


NOTES 


THE CONDUCTIVITY SODIUM LOW TEMPERATURES! 


Woops 


INTRODUCTION 

recent theoretical paper, Ziman (1954) concludes that 
cesses (Peierls 1932; Bardeen 1937) may still play important role the 
resistance the monovalent metals the passage heat and electricity 
quite low temperatures, indeed possibly low ~6/30. suggests that 
higher temperatures the variation with temperature the resistance 
fortuitous agreement with the Bloch theory (which ignores 
cesses) and that only below should the assumptions the Bloch theory 
valid. was therefore decided measure accurately the electrical re- 
sistance sodium very low temperatures because from other evidence 
sodium appears the best approximation ideal monovalent metal, 
has relatively high that the temperature region interest more easily 
accessible, and can obtained extremely pure state. 

Assuming rule, may write the over-all resistivity: 

potpi 


where the residual resistivity due scattering the electrons im- 
purities and the resistivity due thermal scattering. therefore 
small possible. Even with sodium the highest purity available, 
only 1/200 deviations from the Bloch theory, which predicts 
(e.g. Wilson 1953), are detected, then necessary measure with 
error of, most, few parts ten thousand. 

THE EXPERIMENTS 

Two specimens were cast soft glass capillary molds with platinum elec- 
trodes fused in. The capillaries were filled with sodium which was first melted 
vacuum and then forced with small pressure helium gas. The speci- 
mens were later sealed off. Specimen was about cm. long, wound 
The which the measurements were made, together with the method 
temperature control and determination, has been described previously 
White and Woods (1955). 

the measurement the electrical resistance direct reading 
galvanometer amplifier (MacDonald 1947) was used. The temperature the 
specimen was lowered ~4° and current, through stable constant 
resistance, series with the specimen, was adjusted provide voltage 


N.R.C. No. 3839. 
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approximately equal and opposite that across the specimen. This essen- 
tially cancelled the deflection due making possible much more accurate 
determinations p;. The exact fraction the voltage that was due was 
determined later from measurements down Since was stable 
one part 5000 throughout the complete experiment and the current through 
the specimen (100 ma.) was measured least one part 5000, the possible 
error determining should about 5%. This error may in- 
creased perhaps owing inaccuracy reading the small deflections 
(~2 cm.) the galvanometer amplifier. 

the second experiment, gridded galvanometer amplifier (Jones 
1951) was used. The measurements were now limited variations 
the thermal voltages the wires leading fron the cryostat the room 
temperature These were such produce error about 

was found cooling the samples from room temperature below 
that value was observed which remained constant long the tempera- 
ture was not raised subsequently above During one run this period was 
longer than two hours. However the temperature was raised above 10° 
and then returned below the residual resistance was found have in- 
creased and remained the new value long the 
temperature was not raised again. This behavior was observed for each the 
experimental runs. Raising the temperature high 40° between 
readings the residual resistance did not change the effect. Since error 
determining only affected the calculation seriously below about 
K., measurement was made and was then assumed remain 
constant least for the 15-min. period required make the determinations 
between and With this procedure the results all runs give single 
curve for function temperature. 


RESULTS 
Within the experimental error the results may well represented the 
following power laws: 


The experimental values are shown Fig. together with graphs 
the functions given above. Varying the exponents these laws 
much would vitiate the agreement with the experimental data. The 
mean value obtained for the residual resistivity these samples 


DISCUSSION 


From 15° the temperature variation the resistance our speci- 
mens agreement with the Bloch theory and with the results MacDonald 
and Mendelssohn (1948, 1950) and Berman and MacDonald (1951), which 
cover this range, although Berman and MacDonald find absolute values that 
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are somewhat lower. the basis Ziman’s theory, might thought that 
the accord with the Bloch prediction indeed fortuitous, and that the more 
rapid decline the resistance below ~8° due the diminishing 
Umklapp processes the temperature lowered. that case, would 
appear that Umklapp processes are, fact, important temperatures 
low (4° K.) sodium. the other hand the thermal conductivity, 
which was measured the same time low 6/100, shows sign this 
effect. Further discussion this will appear later paper the thermal 
conductivities the alkali metals White, Woods, and MacDonald (1956). 

Recent experiments Barrett (1955) indicate that sodium undergoes 
crystallographic transformation below ~40° The type transformation 
(martensitic) that observed sensitive the stress and purity conditions 
the metal and may exhibit strong hysteresis effects with temperature 
changes. Consequently physical properties measured the temperature region 
such transformation may involve unknown and varying amounts 
different structures and, hence, may not provide reliable test the theory 
these properties. the other hand, the results are interpreted this 
basis, the reproducibility the measurements from run run—especially 
for samples such different physical dimensions—seems surprising. 
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indebted Dr. MacDonald for suggesting this investigation 
and for discussions it. wish thank Dr. Chambers also for his 
interest the problem. 
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STRESS DEFORMATION RELATIONS FOR SOLIDS 
ERICKSEN 


are concerned here with theories solids which the stress assumed 
function displacement gradients, the case numerous theories 
elasticity and plasticity. From the notion that, with some exceptions, 
nonzero amount work required alter the deformation, conclude 
that the stress tensor, tensor proportional it, derivable from scalar 
potential function. This follows easily using theorem Pfaffian equations, 
the proof and first physical applications which are discussed Cara- 
work thermodynamics (1909). Briefly, the theorem that 
Pfaffian equation, 


N 
(1) F (yo, = 0, 


has the property that, every neighborhood arbitrary point there 
exists point such that curve satisfying (1) joins then there 

where are rectangular Cartesian coordinates material particle 
time its initial coordinates. Here and the following, indices take 
values assumption, the stress components are functions the 
derivatives For simplicity, assume there are constraints 
the deformation. 

The rate which the stresses work deforming the material contained 
volume which moves with the material given the well-known 
integral 
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initially occupied the material which occupies time For this rate 
vanish for every volume the differentials must satisfy 
which Pfaffian equation, being from our assumptions function the 
now restrict our attention stress-deformation relations such that, 
arbitrarily close each deformation, there are other deformations which 
unlikely that relations not satisfying this condition can describe real materials. 
Carathéodory’s theorem, applied the Pfaffian considered here, then implies 
the existence functions and such that 


or, equivalently, 


deriving (3) and (4), have not assumed that the processes for which the 
stress-deformation relations hold are reversible that they are either iso- 
thermal adiabatic. The results then apply irreversible well reversible 
processes. When applied Carathéodory’s theorem, tacitly assumed that 
the stress depends only displacement gradients or, one prefers, these 
and other quantities which are held constant while the material deforms. 
The range applicability (3) and (4) should thus not assumed greater 
than the range validity this assumption. 

the special case where and only this case, may write 


where dy. then follows from (2) that the work per unit undeformed 
volume done changing the deformation from the same for all 
motions connecting these, being given This the case 
ordinarily considered nonlinear elasticity, being the so-called elastic 
potential strain energy. discussion this case given Truesdell (1952). 
The existence elastic potential can inferred from thermodynamic 
considerations for reversible isothermal and adiabatic processes. Our result 
stronger the sense that applies wider range phenomena, weaker 
the sense that cannot deduce (5), but only the less-restrictive forms (3) 
and 
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